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ABSTRACT 


It  is  well  recognized  that  conditions  in  the  environment  surround¬ 
ing  ice  nucleu  (IN)  particles  must  be  accurately  controlled  in  order  to 
understand  their  nucleation  behavior.  Concerning  the  condensation¬ 
freezing  and  deposition  mechanisms  of  ice  nucleation,  the  main  factors 
representing  the  environment  are  supersaturation  and  supercooling 
(temperature).  Starting  from  the  concept  of  the  wedge-shaped  ice  ther¬ 
mal  diffusion  chamber  of  Schaller  and  Fukuta,  and  that  of  the  Fukuta- 
Saxena  cloud  condensation  nucleus  (CCN)  spectrometer,  and  by  addition¬ 
ally  considering  operational  efficiency  and  accuracy,  a  continuous-flow, 
horizontal  gradient,  ice  thermal  diffusion  chamber  has  been  developed. 
The  chamber  consists  of  three  parts,  i.e.,  preprocessing,  main  activa¬ 
tion  and  ice  crystal  settling  sections.  A  common  top  plate  is  used  for 
all  three  sections  and  is  held  isothermal  using  a  circulating  bath  while 
a  temperature  gradient  is  maintained  across  the  bottom  plate  of  the  main 
section  by  applying  thermoelectric  modules.  In  the  preprocessing  and 
main  sections,  both  plates  are  covered  with  ice.  A  newly  developed 
method  is  applied  to  smoothly  coat  ice  on  the  plates.  This  configura¬ 
tion  of  the  main  section  results  in  a  nearly  constant  temperature  hori¬ 
zontally  in  the  direction  across  the  sample  flow,  and  produces  a  range 
of  supersaturations.  Heat  pipes  are  utilized  on  the  sides  of  the 
bottom  plate  to  insure  temperature  uniformity  along  the  direction  of 
the  flow.  The  sample  volume  is  sandwiched  in  the  region  of  maximum 
suoersaturation  between  layers  of  filtered  and  predried  air  such  that. 
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as  the  sample  enters  the  main  section,  a  nearly  constant  supersaturation 
is  achieved  vertically  through  the  sample  without  transient  supersatura¬ 
tions.  The  preprocessing  section  is  held  isothermal  at  the  top  plate 
temperature.  The  design  and  flow  of  the  main  section  permit  the  nucle¬ 
ated  ice  crystals  to  be  carried  into  the  ice  crystal  settling  section 
without  loss.  The  bottom  plate  in  the  settling  section  is  maintained  at 
a  temperature  slightly  lower  than  that  of  the  top  plate.  Formation  of 
transient  supersaturations  at  the  entrance  of  the  settling  section  has 
been  avoided  by  delaying  vapor  diffusion  while  allowing  thermal  diffu¬ 
sion  to  proceed.  The  problem  of  transient  supersaturation  development 
has  been  examined  for  shear  (Poiseuille)  flow  cases  using  a  numerical 
method  instead  of  the  commonly  used  slab  flow  assumption.  Wall  effects 
have  also  been  computed  for  the  supersaturation  distribution.  The  sta¬ 
bility  of  the  sample  flow  throughout  the  entire  chamber  has  been  con¬ 
firmed  with  smoke  tests.  A  new  method  for  ice  crystal  detection  has 
been  develoDed  and  is  employed  in  the  ice  crystal  settling  section. 

Mylar  copy  film  (carbon  paper)  holding  condensed  water  droplets  is 
placed  in  the  settling  section.  Ice  crystals  nucleated  in  the  main 
section  fall  on  the  film  and  grow  to  visible  sizes  in  the  presence  of 
the  droplets.  The  positions  ofthe  ice  crystals  across  the  flow  direc¬ 
tion  give  the  supersaturations  and  temperatures  at  which  they  nucleated. 
The  optimum  design  of  the  chamber  in  terms  of  the  height,  width,  length 
and  flow  speed  has  been  determined. 

Two  numerical  models  have  been  developed  to  compute  the  temperature 
and  supersaturation  fields  within  the  chamber.  The  first  computes  the 
steady  state  conditions  in  a  vertical  plane  perpendicular  to  the  flow 


and  includes  the  wall  effects.  The  second  computes  the  steady  state 
conditions  in  a  vertical  plane  parallel  to  the  flow  using  the  shear  flow 
velocity  profile  and  the  total  pressure-  and  temperature-dependent  co¬ 
efficients  for  vapor  and  thermal  diffusion. 

The  increased  accuracy  of  the  developed  chamber  has  resulted  in 
several  significant  new  findings.  For  all  the  temperature  and  supersa¬ 
turation  ranges  studied,  the  rate  of  deposition  nucleation  is  much 
lower  than  that  for  condensation-freezing.  For  silver  iodide  (Agl) 
particles  of  average  diameter  0.05  urn,  only  a  small  fraction  are  active 
below  1%  water  supersaturation,  while  for  1,5  -  dihydroxynaphthalene 
(DN)  of  average  diameter  0.1  um,  a  much  larger  fraction  are  active. 

The  total  number  of  active  nuclei  for  these  samples  was  approximately 
three  times  greater  for  ON  than  for  Agl  at  high  supersaturations.  For 
temperatures  warmer  than  -10°C,  the  condensation-freezing  mechanism  is 
not  effective  for  Agl.  Natural  nuclei  and  kaolinite  show  little  super¬ 
saturation  dependence  below  1%  water  supersaturation  and  as  much  as  an 
order  of  magnitude  increase  as  water  supersaturation  increases  from  1% 
to  2.5%.  Additionally  for  natural  nuclei,  large  variances  were  observed 
with  high  counts  in  cloud  downdrafts  and  low  counts  during  precipitation. 
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CHAPTER  1 


INTRODUCTION 

The  occurrance  of  precipitation  in  the  atmosphere  has  a  large  im¬ 
pact  on  man's  activities,  often  dictating  where  civilizations  locate, 
how  they  evolve  and  even  whether  they  survive.  The  absence  of  precipi¬ 
tation  brings  famine  and  despair  whereas  an  over-abundance  can  also 
lead  to  similiar  results.  Mankind  has  to  a  large  extent  adapted  to  the 
global  distribution  of  precipitation  by  migrating  to  regions  where  pre¬ 
cipitation  is  adequate  and  avoiding  those  where  moisture  supplies  are 
inadequate.  However,  even  within  regions  of  normally  adequate  precipi¬ 
tation,  short  or  long  term  fluctuations  can  result  in  drought  or  flood¬ 
ing  and  thereby  exert  significant  effects  on  man’s  environment.  Both 
the  geographic  and  time  variances  are  to  a  large  extent  a  result  of 
fluctuating  moisture  availability  in  the  environment.  However,  the 
mechanisms  involved  in  the  process  of  converting  the  water  suspending 
in  the  atmosphere  into  precipitation  elements  also  play  significant 
roles.  Hence,  the  study  of  these  mechanisms  involved  in  precipitation 
formation  is  of  great  importance  to  man. 

Precipitation  transports  moisture  from  the  atmosphere  to  the 
earth's  surface  in  various  forms  including  rain,  snow,  hail  and  graupel 
(Peterssen,  1958).  It  is  important  to  note  that  with  the  exception  of 
some  tropical  showers,  the  ice  phase  is  always  involved  in  the  precipi¬ 
tation  formation  process  even  when  the  final  form  observed  on  the 
ground  is  rain  (Mason,  1971).  These  raindrops  form  when  ice  crystals 
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grow  in  supercooled  clouds  (i.e.  clouds  consisting  of  water  droplets 
below  0°C)  and  subsequently  fall  and  melt  at  lower  altitudes.  This 
mechanism  of  precipitation  formation,  known  as  the  Bergeron  mechanism, 
plays  a  dominant  role  in  the  atmosphere  (Bergeron,  1935).  It  is  the 
process  of  the  initial  ice  crystal  formation  involved  in  this  mechanism 
that  receives  primary  attention  in  this  study. 

Formation  of  ice  crystals  or  ice  nucleation  must  take  place  at 
least  at  the  early  stage  and  possibly  throughout  the  Bergeron  process 
of  precipitation  formation.  Ice  crystal  nucleation  can  be  categorized 
into  two  types,  i.e.  homogeneous  and  heterogeneous.  In  the  former,  the 
phase  change  to  the  solid  (crystal)  state  occurs  without  the  assistance 
of  a  foreign  substance  while  the  latter  involves  some  foreign  particle 
or  substrate  in  the  phase  change  process.  Such  particles  are  called 
ice  nuclei.  Since  temperatures  of  -40°C  or  lower  are  required  for 
homogeneous  ice  nucleation  (Shaefer,  1948),  the  occurrence  of  signifi¬ 
cant  numbers  of  natural  ice  crystals  in  atmospheric  processes  at  much 
warmer  temperatures  indicates  that  this  is  not  the  primary  mechanism  in 
nature  in  most  cases  and  in  fact  formation  of  the  cirrus  anvils  asso¬ 
ciated  with  thunder  or  hail  clouds  is  probably  one  of  the  rare  examples 
of  the  process  in  the  atmosphere.  Th.is  leaves  heterogeneous  nucleation 
as  the  contributing  mechanism  for  ice  crystal  formation  for  most  cases 
in  nature  (Weickmann,  1971;  Rodgers,  1976). 

Investigators  studying  heterogeneous  ice  nucleation  in  the  atmos¬ 
phere  have  determined  that  there  are  large  fluctuations,  both  in  time 
and  space,  in  the  number  of  particles  which  act  as  ice  nuclei  and  fur¬ 
thermore  that  their  number  is  normally  quite  small  compared  to  the 
number  of  cloud  condensation  nuclei  (i.e.  nuclei  which  aid  the  formation 
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of  cloud  droplets)(Fukuta,  1966;  Weickmann,  1971).  Three  orders  of 
magnitude  differences  are  quite  common  between  the  ice  and  the  cloud 
condensation  nuclei.  This  scarcity  of  natural  ice  nuclei  leads  investi¬ 
gators  to  the  realization  that  opportunities  may  exist  to  enhance  the 
precipitation  process  by  artificially  adding  effective  ice  nuclei  to  a 
supercooled  cloud  environment.  Therefore,  the  study  of  ice  nucleation 
in  the  atmosphere  concentrates  on  the  effective  number  and  activation 
characteristics  of  heterogeneous  ice  nuclei  both  naturally  occurring 
and  artificially  generated.  For  ice  nuclei  existing  naturally  in  the 
atmosphere,  the  studies  have  had  very  limited  success  due  mainly  to  the 
scarcity  and  variability  of  these  particles.  Cloud  physicists  have 
found  materials  which,  when  dispersed  properly  into  a  supercooled  cloud 
environment,  behave  as  effective  heterogeneous  ice  nuclei.  These  in¬ 
clude  silver  iodide,  lead  iodide,  metaldehyde  and  1  ,5-dihydroxynaphtha- 
lene.  These  ice  nuclei  have  been  studied  extensively  but  findings  have 
been  inconsistent  due  primarily  to  the  differences  in  the  methods  of 
measurement  reflected  in  the  designs  of  instrumentation  used  (Hallett, 
1971;  Schaefer,  1971;  Langer,  1973;  Vali,  1976). 

The  previous  studies  have  indicated  that  there  are  three  distinct 
mechanisms  for  heterogeneous  ice  nucleation  (Schaller  and  Fukuta,  1979). 
The  first,  deposition  nucleation  involves  the  direct  transition  from 
the  vapor  to  the  solid  state  with  some  foreign  particles  acting  as  the 
nuclei  for  crystallization.  The  condensation-freezing  nucleation 
mechanism  involves  all  three  states,  with  vapor  condensing  around  the 
nucleus  to  form  a  liquid  droplet  which  supsequently  freezes.  The  final 
mechanism,  contact-freezing  nucleation,  also  produces  a  crystal  in  the 
supercooled  droplet  but  the  nucleus  is  external  to  the  droplet  and 


initiates  freezing  upon  contact  with  the  droplet  surface.  Each  mechan¬ 
ism  has  its  own  dependencies  on  temperature,  supersaturation,  and  time 
as  well  as  on  the  characteristics  of  the  nuclei  involved. 

Researchers  have  used  many  different  methods  in  their  study  of 
ice  nuclei  but  generally  they  can  be  grouped  into  three  categories, 
the  cloud  chamber  method,  the  filter  method  and  the  large  crystal  method. 
Most  mixing  chambers  sustain  supercooled  clouds  in  them  as  a  reproduc¬ 
tion  of  atmospheric  conditions  and  permit  studying  the  behavior  of  ice 
nuclei  injected  in  the  clouds.  There  are  a  few  problems  associated 
with  this  method.  First,  during  the  processes  of  sample  ice  nuclei  in¬ 
troduction  and  humidification,  transient  supersaturation  is  likely  to 
occur,  making  the  condition  of  ice  nuclei  activation  unclear.  Second, 
only  one  supersaturation  level,  i.e.  water  saturation,  can  be  stably 
sustained.  Third,  ice  nuclei  behaviors  under  different  levels  of  super¬ 
saturation  cannot  be  studied.  Different  levels  of  supersaturations  can 
be  produced  by  expansion  chambers  but  they  cannot  be  sustained  under 
steady  state.  In  addition,  ice  crystal  counts  obtained  are  results  of 
different  ice  nucleation  mechanisms  occurring  simultaneously. 

The  filter  method  allows  the  ice  nuclei  removed  from  the  air  to  be 
activated  under  different  supersaturations  and  temperatures.  Results 
obtained  by  this  method,  however,  suffer  from  effects  of  coexisting 
impurities  and  substrate. 

While  the  clean  faces  of  large  cleaved  single  crystals  of  an  ice 
nucleus  compound  furnish  excellent  macroscopic  surfaces  for  studying 
ice  nucleation  mechanisms,  the  important  size  effect  of  the  nuclei  is 
lost  in  addition  to  unrealistically  magnified  effect  of  steps  and  cavi¬ 
ties  (Schaller  and  Fukuta,  1979). 
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Although  these  methods  have  different  advantages  and  disadvantages 
in  heterogeneous  ice  nucleation  studies,  the  cloud  chamber  method  can 
avoid  the  substrate  problem  and  does  include  the  size  effect.  There¬ 
fore,  use  of  a  cloud  chamber  would  give  the  best  opportunity  to  gather 
reliable  and  consistent  data,  had  the  difficulty  of  changing  and  main¬ 
taining  the  supersaturation  at  different  levels  been  solved. 

Thermal  diffusion  cloud  chambers  are  capable  of  stably  maintaining 
temperature  and  supersaturation  at  any  combined  levels,  thereby  elimi¬ 
nating  the  above  mentioned  deficiency  of  cloud  chambers.  The  existing 
cloud  chambers  of  this  type,  however,  can  only  accommodate  small  sample 
volumes,  hold  one  set  of  temperature  and  supersaturation  at  the  desired 
portion  of  the  chambers,  and  the  counting  schemes  employed  are  unreli¬ 
able  (normally  visual  method).  The  wedge-shaped  ice  thermal  diffusion 
chamber  recently  introduced  by  Schaller  and  Fukuta  (1979)  produces  a 
range  of  supersaturation  and  therefore  eases  the  disadvantage  of  single 
set  of  temperature  and  supersaturation  in  the  normal  thermal  diffusion 
cloud  chambers.  However,  it  still  suffers  from  the  small  sampling  vol¬ 
ume.  In  addition,  the  time  dependency  of  ice  nucleation  cannot  be 
accurately  studied.  The  purpose  of  this  study  is  to  design  and  con¬ 
struct  a  new  horizontal  gradient,  continuous  flow,  ice  thermal  diffu¬ 
sion  cloud  chamber  which  maintains  the  assets  of  the  present  chambers 
while  eliminating  some  of  the  difficulties.  This  new  cloud  chamber  is 
then  used  to  study  various  ice  nuclei  and  identify  their  nucleation 
mechanisms  as  well  as  their  dependence  on  temperature  and  supersatura¬ 
tion. 

A  review  of  the  previous  methods  of  ice  nuclei  study  is  presented 
in  Chapter  2.  Chapter  3  includes  the  design  and  construction  of  the 
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new  chamber  with  a  discussion  of  the  distinct  improvements  over  pre¬ 
vious  chambers.  The  numerical  models  used  to  simulate  the  conditions 
in  the  chamber  are  included  in  Chapter  4.  Chapter  5  presents  the  pro¬ 
cedure  of  chamber  operation  and  sampling.  Results  and  discussions  are 
included  in  Chapter  6  with  conclusions  in  Chapter  7. 


CHAPTER  2 


METHODS  USED  IN  PREVIOUS  ICE  NUCLEATION  STUDIES 

Since  the  primary  purpose  for  studying  ice  nuclei  is  to  understand 
their  behaviors  under  naturally  existing  conditions,  it  is  highly  de¬ 
sirable  to  reproduce  the  conditions  as  faithfully  as  possible 
(Podzimek,  1971).  For  studies  performed  under  unnatural  conditions, 
the  results  must  be  manipulated  to  compensate  for  these  deviations. 

Any  instrument  which  can  accurately  control  the  environment  within  the 
restraints  of  nature  will  therefore  have  a  significant  advantage  over 
instruments  where  unnatural  conditions  exist,  even  if  only  for  short 
periods  of  time. 

For  ice  nucleation  studies,  temperature  has  always  been  an  impor¬ 
tant  parameter  and  has  been  fairly  accurately  controlled  in  most  stud¬ 
ies  (Mason  and  Hallett,  1956;  Bigg,  1957).  The  more  recent  recognition 
of  the  roles  of  supersaturation  in  heterogeneous  ice  nucleation  has 
clarified  the  importance  of  its  accurate,  though  difficult,  control 
(Hudson  and  Squires,  1976;  Schaller  and  Fukuta,  1979).  The  techniques 
used  in  ice  nucleation  studies  are  surveyed  below  and  should  be  criti¬ 
cally  evaluated  for  their  control  of  temperature  and  supersaturation  in 
order  to  determine  the  usefulness  of  the  results  when  applied  to 
realistic  atmospheric  conditions. 

2. 1  Cloud  Chambers 

Although  there  are  a  variety  of  types  of  cloud  chambers,  their  basic 
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purpose  is  to  furnish  an  enviroment  of  a  supercooled  cloud  in  which 
ice  nucleation  can  be  studied.  Schaefer  (1946)  first  used  a  small  com¬ 
mercial  freezer  as  a  cloud  chamber  by  blowing  breath  into  it,  and  since 
then  other  investigators  have  developed  various  types  of  cold  chambers, 
each  with  its  own  distinct  advantages  and/or  disadvantages. 

2.1.1  Mixing  chambers 

The  walls  of  mixing  chambers  are  normally  maintained  at  a  constant 
temperature  with  the  supercooled  cloud  formed  either  by  cooling  the  sam¬ 
ple  air  below  its  dew  point  temperature  or  by  the  addition  of  water  va¬ 
por  or  droplets.  The  intent  is  to  provide  an  isothermal  environment 
with  a  constant  supersaturation  field;  but  in  reality,  cloud  tempera¬ 
tures  can  vary  up  to  5°C  (Bigg,  1957)  and  large  fluctuations  in  super¬ 
saturation  often  occur  (Ohtake,  1976).  In  addition,  selection  of 
steady  state  supersaturation  at  an  arbitrary  level  is  not  possible,  and 
the  values  which  do  occur  cannot  be  accurately  determined.  There  also 
exists  the  possibility  of  transient  supersaturations  (i.e.  supersatura¬ 
tions  of  short  duration  which  are  much  greater  than  the  final  values). 
These  tend  to  result  in  unrealistically  high  numbers  of  active  nuclei. 

2.1.2  Expansion  chambers 

Expansion  chambers  were  originally  used  in  nuclear  physics  to 
trace  cosmic  rays  (Sears  and  Zemansky,  1962)  and  later  adapted  for  use 
in  cloud  physics.  For  ice  nucleation  studies,  the  air  is  normally 
cooled  to  a  subfreezing  temperature  and  then  further  cooled  by  rapid 
adiabatic  expansion  (Warner,  1957).  The  results  are  reasonably  repro- 
duciable  for  a  given  instrument,  but  a  significant  variance  has  been 


noted  between  different  chambers,  even  ones  with  identical  design 
(Kline  and  Brier,  1961).  Expansion  chambers  are  reliable  and  readily 
adaptable  for  both  field  and  laboratory  use  but  suffer  from  some  basic 
disadvantages.  The  expansion  produces  short  term  supercooling  and  su¬ 
persaturation,  which  can  be  estimated  with  reasonable  accuracy,  but  the 
effect  of  the  short-lived  conditions  cannot  be  exactly  evaluated.  Nor¬ 
mally,  the  transient  supersaturations  produced  in  the  chamber  are  much 
higher  than  those  which  occur  in  nature,  resulting  in  unreliably  high 
nuclei  counts. 

2.1.3  Diffusion  chambers 

Langsdorf  developed  a  diffusion  cloud  chamber  for  use  in  nuclear 
physics  in  1936.  He  suggested  that  it  might  be  useful  for  atmospheric 
studies,  but  it  was  not  utilized  for  such  work  until  Schaefer  developed 
several  types  for  atmospheric  use  (Schaefer,  1952,  1954).  These  cham¬ 
bers  are  convectively  stable  and  are  capable  of  reproducing  tempera¬ 
tures  and  supersaturation  fields  with  excellent  accuracy  (Squires, 
1972).  Caution  must  be  taken  to  avoid  wall  effects  and  transient 
supersaturations  while  ensuring  that  the  counting  method  employed 
yields  accurate  results.  On  the  other  hand,  thermal  diffusion  chambers 
suffer  from  the  disadvantage  of  necessarily  having  small  heights  and 
consequently  small  volumes  (Squires,  1971),  thereby  limiting  their  use 
for  samples  with  a  small  ice  nuclei  population  (Hallett,  1971). 

2.2  Precipitation  Methods 

Instruments  using  this  method  collect  the  ice  nuclei  on  substrate; 
and  subsequently  activate  them  under  controlled  temperatures  and 
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supersaturations  and  then  grow  the  crystals  to  detectable  size.  Al¬ 
though  these  methods  yield  information  about  the  sizes  of  the  nuclei, 
they  are  critically  dependent  on  the  development  procedure  and  suffer 
from  a  lack  of  reproducibility  (Langer,  1973).  Furthermore,  the  effect 
of  the  substrate  can  contribute  significantly  to  the  apparent  activa¬ 
tion  characteristics  of  the  nuclei. 

2.2.1  Impaction  method 

Ice  nuclei  particles  are  allowed  to  impact  on  some  surface  and  are 
subsequently  developed.  Impaction  on  a  varnished  surface  was  used  by 
Yang  (1962)  while  Gerber  (1976)  used  a  foil  in  the  Goetz  aerosol  spec¬ 
trometer,  The  collection  efficiency  of  these  methods  is  somewhat  un¬ 
certain  and  the  substrate  problem  is  certainly  present. 

2.2.2  Membrance  filter  method 

This  method  makes  possible  the  sampling  of  large  volumes  of  air 
with  a  relatively  simple  experimental  setup.  The  method  of  using 
millipore  filters  was  first  employed  by  Bigg,  et  al . ,  (1961)  and  has 
subsequently  been  used  by  many  investigators  (Bigg,  et  al . ,  1963).  The 
method,  however,  is  dependent  on  the  type  of  filter  used,  the  amount  of 
air  drawn  through  it,  and  the  development  process.  Again,  significant 
substrate  effect  is  expected  with  this  method. 

2.3  Use  of  Thermal  Diffusion  Chambers  for  Cloud  Condensation 

Nuclei  (CCN)  Studies 

Of  all  the  types  of  instruments  used  in  atmospheric  nuclei  studies, 
the  thermal  diffusion  chamber  is  unique  in  its  ability  to  produce  ac¬ 
curately  controlled  temperature  and  supersaturation  fields  under  steady 
state  while  avoiding  the  troublesome  substrate  effects  and  incoporating 
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the  size  effect.  These  chambers  operate  with  a  thermally  stable  inter¬ 
ior  environment  where  the  temperature  profiles  can  be  precisely  deter¬ 
mined.  As  is  the  case  in  all  other  ice  nucleus  detection  methods,  the 
supersaturation  values  cannot  be  directly  measured.  However,  with 
known  boundary  conditions  for  vapor  pressure  and  temperature,  the  super 
saturation  fields  can  be  accurately  computed.  Adequate  knowledge  of 
these  parameters  permits  detailed  studies  of  ice  nuclei  to  be  performed 
The  thermal  diffusion  chambers  used  in  atmospheric  studies  have 
taken  many  different  designs.  Concentric  cylinders  with  wetted  wall 
have  been  used  (Rosen  and  Hofmann,  1977)  although  the  most  common  de¬ 
sign  utilizes  parallel  plates  (Katz  and  Ostermier,  1967;  Radke  and 
Hobbs,  1969;  Elliott,  1971;  Heist  and  Reiss,  1973).  While  most  of  the 
thermal  diffusion  chambers  process  the  sample  under  static  conditions, 
in  recent  years,  continuous  flow  thermal  diffusion  chambers  have  been 
developed.  Table  1  lists  the  use  of  thermal  diffusion  chambers  in  the 
study  of  atmospheric  CCN  (Schaller,  1975). 

2.4  Use  of  Thermal  Diffusion  Chambers  for  Ice  Nuclei  (IN)  Studies 

By  lowering  the  temperature  and  substituting  ice  surfaces  for  the 
wetted  ones,  the  thermal  diffusion  chambers  used  for  CCN  studies  can  in 
principle  be  applied  to  ice  nuclei  studies.  Diffusion  chambers  were 
used  in  IN  studies  to  develope  filters  as  discussed  earlier  in  addition 
to  investigating  ice  nucleation  both  on  large  crystal  surfaces  and  on 
particles  freely  suspended  in  air. 

2.4.1  For  studies  using  large  crystal  surfaces 

By  placing  a  crystal  surface  in  the  interior  of  a  thermal  diffu¬ 
sion  chamber,  the  environment  over  the  surface  can  be  controlled 
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precisely  and  ice  nucleation  can  be  studied  in  detail.  Unfortunately , 
the  results  of  these  studies  cannot  be  applied  directly  to  atmospheric 
processes  since  the  size  effect  is  not  incorporated  (Fletcher,  1958) 
and  the  effect  of  the  large  substrate  is  not  known. 

2.4.2  For  studies  of  freely  suspended  nuclei 

The  attributes  of  the  thermal  diffusion  chamber  method  have  been 
utilized  very  effectively  by  Schaller  and  Fukuta  (1979)  in  a  wedge- 
shaped  chamber.  The  isothermal  plates  were  replaced  with  ones  where  a 
temperature  gradient  is  maintained  across  the  plates,  thereby  producing 
a  center  environment  where  a  range  of  supersaturations  exist.  This 
chamber  allowed  the  detailed  study  of  artificial  ice-nucleating  agents, 
but  the  limited  volume  precluded  its  use  for  natural  aerosols.  The 
visual  counting  method  leaves  some  room  for  error,  and  wall  effects  may 
affect  the  interior  environment. 

2.5  Principle  of  Operation  of  Thermal  Diffusion  Chamber 

Since  the  thermal  diffusion  chamber  method  is  capable  of  the  accu¬ 
rate  temperature  and  supersaturation  control  while  avoiding  substrate 
effects  and  including  size  effects,  it  has  the  best  promise  for  use  in 
detailed  ice  nuclei  studies.  The  mechanism  whereby  ice  and  water  super¬ 
saturations  are  sustained  in  subfreezing  temperature  ranges  is  illus¬ 
trated  in  Figure  1.  The  thermal  diffusion  chamber  works  by  evaporating 
water  vapor  from  a  saturated  warm  surface  (the  top  plate),  diffusing  it 
through  the  chamber  interior  and  condensing  it  on  a  cooled  surface  (the 
bottom  plate).  The  top  plate  is  ice  covered  and  maintained  at  a  tem¬ 
perature  T^.  The  bottom  plate  is  also  ice  ccvered  and  maintained  at 
another  temperature  T^  where  T^T^..  Katz  and  Mirabel  (1  975)  showed 


VAPOR  PRESSURE  (mb) 


Figure  1.  Principle  of  supersaturation  development  in  thermal 
diffusion  chambers.  Point  A  is  at  the  top  plate 
temperature,  C  at  the  bottom  plate  temperature,  with 
point  B  representing  the  conditions  at  the  midpoint 
between  the  plates. 
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that  the  temperature  and  vapor  pressure  profiles  between  the  plates 
are  very  close  to  linear.  This  profile  is  shown  in  Figure  1  by  the 
straight  line  from  A  to  C.  However,  the  saturation  values  for  both  ice 
and  water  form  curves  which  are  concave  upward  with  the  curve  for  water 
saturation  slightly  above  that  for  ice.  Hence,  for  ice  coated  bounda¬ 
ries  where  the  temperature  is  known,  the  vapor  pressure  can  be  accurate¬ 
ly  determined  using  the  ice  saturation  values.  In  Figure  1,  the  top  and 
bottom  plate  values  are  given  by  points  A  and  C,  respectively.  The 
linear  profile  through  the  interior  results  in  vapor  pressure,  P,  higher 
than  ice  saturation,  P  and  ice  supersaturation  results.  Ice  super¬ 
saturation  is  defined  as  S.  =  (P/P  .)-l .  For  sufficiently  large  temper¬ 
ature  differences,  vapor  pressure  values  greater  than  water  saturation, 

Pc  .  exist  and  water  supersaturated  conditions  are  created.  Water  su- 
sw 

persaturation  is  defined  as  Sw  =  (P/Psw)-1-  These  temperature  and  va¬ 
por  pressure  profiles  produce  maxima  for  both  water  and  ice  supersatur¬ 
ation  around  the  chamber  center,  but  they  are  not  necessarily  colocated. 

Most  thermal  diffusion  chambers  in  use  have  isothermal  plates  which 
produce  the  same  temperature  and  vapor  pressure  profiles  throughout  the 
chamber  interior  except  for  the  regions  close  to  the  walls.  The  chamber 
center  is  maintained  at  a  constant  temperature  and  supersaturation.  The 
CCN  spectrometer  of  Fukuta  and  Saxena  (1979a)  as  well  as  the  wedge- 
shaped  ice  thermal  diffusion  chamber  of  Schaller  and  Fukuta  (1979)  in¬ 
corporate  temperature  gradients  on  the  plates,  thereby  creating  a  range 
of  supersaturations  horizontally  across  the  chamber  center.  This 
allows  a  sample  to  be  studied  under  a  range  of  supersaturations  during 
an  experimental  run  instead  of  under  only  a  single  value  of  supersatur¬ 
ation.  This  not  only  expedites  the  sampling  procedure  but  minimizes 
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the  errors  caused  by  changes  in  the  sample  between  experiments. 

The  ice  thermal  diffusion  chamber  presented  in  this  study  is  some¬ 
what  different  from  the  two  mentioned  above  in  that  the  top  plate  is 
held  isothermal  while  a  temperature  gradient  is  maintained  across  the 
bottom  plate.  The  result  is  again  a  range  of  supersaturation  across 
the  chamber  center  with  an  accompanying  small  temperature  change.  Fig¬ 
ure  2  illustrates  a  typical  case  where  the  two  lines  ABC  and  ADE  in 
part  (a)  show  the  extremes  which  occur  in  the  sample.  The  positioning 
of  the  sample  will  be  explained  in  the  next  chapter.  It  should  be 
noted  in  Figure  2(b)  that  water  supersaturation  exists  only  in  a  por¬ 
tion  of  the  chamber  while  ice  supersaturation  exists  throughout  the 
interior. 
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CHAPTER  3 


APPARATUS 

3.1  Factors  to  be  Considered  in  the  Design 

Any  instrument  used  for  the  detailed  study  of  ic£  nuclei  behavior 
must  be  capable  of  accurately  controlling  both  the  temperature  and  the 
supersaturation  of  the  environment  in  which  the  nuclei  are  activated 
since  ice  nucleation  is  normally  a  function  of  both  variables  (Schaller 
and  Fukuta,  1979).  As  discussed  earlier,  thermal  diffusion  chambers 
allow  accurate  control  of  these  variables  if  the  boundaries  inside  the 
chambers  can  be  stably  controlled.  Temperature  control  of  the  top  and 
bottom  plates  of  the  chambers  can  be  achieved  by  various  means  includ¬ 
ing  circulating  refrigerated  baths,  thermoelectric  modules  and  heating 
wires.  Accurate  temperature  measurement  is  equally  important  and  can 
be  achieved  with  the  use  of  thermocouples  and  high  accuracy  meters. 

The  high  accuracy  required  for  vapor  pressure  at  the  surfaces  of  top 
and  bottom  plates  may  be  achieved  with  the  use  of  ice  covered  surfaces 
under  controlled  temperatures. 

A  diffusion  chamber  with  isothermal  plates  provides  a  central 
environment  with  only  a  single  set  of  temperature  and  supersaturation. 
By  producing  a  temperature  gradient  on  one  or  both  plates,  a  range  of 
supersaturations  and/or  temperatures  can  be  produced  (Schaller  and 
Fukuta,  1979;  Fukuta  and  Saxena,  1979a  and  1979b).  Simultaneous  expo¬ 
sure  of  ice  nuclei  to  this  range  of  supersaturation  and/or  temperature 
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makes  possible  a  much  more  efficient  and  comprehensive  study  cf  the  ice 
nuclei.  In  this  manner,  ice  nuclei  behaviors  on  a  line  instead  of  a 
single  point  on  the  supersaturation-temperature  plane  can  be  determined 
from  a  single  sampling. 

Most  thermal  diffusion  chambers  have  a  small  height  and  conse¬ 
quently  a  small  volume.  Ice  nuclei  studies  which  require  large  smaple 
volumes  cannot  utilize  such  chambers.  However,  by  creating  a  flow 
through  the  chamber,  studying  large  volume  of  sample  becomes  possible 
(Sinclair  and  Hoopes,  1975).  Caution  must  be  taken  to  insure  the  flow 
is  laminar  so  that  the  environment  where  ice  nuclei  activate  can  be  ac¬ 
curately  determined.  Any  turbulence  in  the  flow  would  mix  the  sample 
air  and  the  clarity  of  the  condition  for  sample  activation  would  be  lost. 

Wall  effects  (Elliott,  1971;  Goroch  and  Carstens,  1972;  Tomlinson 
and  Fukuta,  1979)  and  transient  supersaturations  (Fitzgerald,  1970,1972, 
Saxena  et  al.,  1970)  are  important  in  designing  a  diffusion  chamber  with 
well-defined  chamber  conditions.  Wall  effects  can  be  avoided  by  main¬ 
taining  a  large  aspect  ratio  (width  to  height  ratio)  and  avoiding  the 
region  close  to  the  wall.  Transient  supersaturations  must  likewise  be 
avoided  since  the  momentarily  high  supersaturations  activate  some  ice 
nuclei  which  otherwise  would  remain  unactivated.  Preconditioning  of 
the  sample  can  eliminate  transient  supersaturations  in  many  cases 
(Saxena,  et  al.,  1970).  However,  preconditioning  is  not  always  poss¬ 
ible  and  a  delay  of  the  vapor  diffusion  can  be  used  instead  to  avoid 
transient  supersaturations  (Fukuta  and  Saxena,  1979a). 

Once  the  ice  nuclei  produce  ice  crystals,  care  must  be  taken  that 
all  the  nucleated  crystals  are  counted.  Ice  crystal  counting  in  thermal 
diffusion  chambers  has  been  one  of  the  crucial  problems  to  be  solved. 
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(Hallett,  1971).  The  new  method  developed  in  this  study  allows  the 
crystals  to  be  integrated  on  a  special  film.  It  is  imperative  that  all 
the  nucleated  crystals  be  collected  and  not  allowed  to  melt  or  sublimate 
prior  to  counting.  It  is  equally  important  that  no  crystals  will  be 
produced  by  the  collecting  surface,  i.e.  the  film  will  produce  no  back¬ 
ground  crystals.  By  taking  these  precautions,  we  insure  that  all  the 
crystals  seen  on  the  collection  film  were  nucleated  in  the  chamber 
under  the  given  environmental  conditions. 

Control  and  measurement  of  the  air  flow  is  another  critical  factor. 
Since  the  quantity  to  be  determined  by  measurements  is  the  number  of 
active  ice  nuclei  per  unit  volume,  the  volume  of  the  sample  must  be 
accurately  determined.  The  residence  time  of  the  sample  also  needs  to 
be  computed  using  the  flow  rate,  because  it  is  this  residence  time  that 
dictates  where  and  when  steady  state  is  reached  in  the  chamber. 

The  new  continuous  flow  horizontal  gradient  ice  thermal  diffusion 
chamber  presented  here  satisfies  all  the  requirements  stated  above. 

The  chamber  is  designed  with  three  primary  sections,  i.e.,  the  precham¬ 
ber,  the  main  chamber,  and  the  postchamber  together  with  an  intake  de¬ 
vice  for  proper  introduction  of  the  sample  and  an  exit  device  for  re¬ 
moving  the  air  from  the  chamber. 

3.2  Final  Design  and  Construction 

The  chamber  was  constructed  entirely  in  our  laboratory.  Copper 
was  chosen  as  the  basic  metal  material  for  top  and  bottom  plates  of  the 
chamber  for  its  good  heat  conduction,  ease  of  fabrication  by  soldering, 
and  rustproof  properties.  Plexiglass  was  utilized  where  an  insulating 
material  was  required.  It  was  selected  for  its  ease  of  fabrication, 
reasonable  strength,  and  transparency  where  observation  into  the 
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chamber  was  desirable.  Figure  3  shows  the  basic  design.  Photographs  of 
the  entire  chamber  and  support  equipment  are  presented  in  Figures  4(a) 
and  4(b).  The  main  parts  are  identified  in  Figure  5.  In  the  following 
discussion,  directions  (left  and  right)  are  defined  as  facing  downstream 
of  the  flow. 

3.2.1  Sample  intake  device 

The  sample  should  be  introduced  and  moved  through  a  region  of 
approximately  constant  temperature  and  supersaturation  near  the  vertical 
center  of  the  chamber  while  avoiding  the  regions  close  to  the  walls 
(Hudson  and  Squires,  1973).  This  positioning  is  achieved  by  sandwich¬ 
ing  the  sample  with  filtered  air,  i.e.  air  filtered  of  both  cloud  con¬ 
densation  and  ice  nuclei.  The  filtered  air  also  separates  the  sample 
from  the  walls.  This  sandwiching  of  the  sample  is  done  with  an  intake 
device  which  is  placed  at  the  upstream  end  of  the  chamber.  A  photo¬ 
graph  is  presented  in  Figure  6(a). 

The  filtered  air  and  the  sample  air  enter  tne  intake  device  from 
both  sides  to  insure  the  horizontal  uniformity  of  the  layers  as  shown 
in  Figure  6(b).  Low  angle  wedges  (fifteen  degrees  or  less)  are  used 
for  filtered  air  introduction  into  the  chamber  to  avoid  turbulence. 

Brass  plates  of  1/32  inch  (0.8  mm)  thickness  are  used  in  tne  construc¬ 
tion  and  are  glued  with  epoxy  resin  instead  of  soldered  to  avoid  warp- 
age.  Shims  of  the  appropriate  thickness  were  used  to  obtain  uniformity 
of  the  slit  width  during  the  construction.  The  slit  openings  are  1/16 
inch  (1.6  mm)  for  the  filtered  air  and  1/32  inch  (0.8  mm)  for  the  sam¬ 
ple  air.  The  sides  are  made  of  1/8  inch  (3.2  mm)  plexiglass,  the  right 
side  keot  clear  for  observation  and  the  left  side  painted  with  flat 
fclack.  The  vertical  position  of  the  sample  air  can  be  controlled  by 
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Figure  4(a).  Photograph  of  the  entire  chamber  with  the  top  removed 
and  inverted  behind  the  chamber. 


Figure  4(b).  Photograph  of  the  support  equipment  including  the  flow¬ 
meters,  the  power  supplies,  the  air  pump,  and  the  Terak 
computer. 
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Figure  6(a).  Photograph  of  the  intake  device  before  assembly  showing 
the  wedges  used  for  sandwiching  the  sample  air  between 
the  layers  of  filtered  air. 

Figure  6(b).  Photograph  of  the  intake  device  before  assembly  showing 
the  symmetry  used  for  the  introduction  of  the  filtered 
air  and  the  sample  to  insure  uniformity  of  the  layers. 

Figure  6(c).  Photograph  of  a  smoke  test  showing  the  uniformity  and 
laminar  flow  of  the  sample  air  as  it  is  sandwiched 
between  layers  of  clear  filtered  air. 

Figure  6(d).  Photograph  of  a  smoke  test  viewed  from  the  exit  end  of 
the  postchamber. 
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adjusting  the  ratio  of  the  thicknesses  of  the  filtered  air  above  and 
below  the  sample.  The  thickness  of  the  sample  layer  is  adjusted  by 
controlling  its  flow  rate.  The  sample  layer  is  separated  from  the  walls 
by  blocking  the  slit  for  approximately  2.5  cm  adjacent  to  each  wall. 
Smoke  was  used  to  verify  the  position  and  uniformity  of  the  sample  layer. 
Photographs  of  smoke  tests  are  shown  in  Figures  6(c)  and  6(d). 

3.2.2  Prechamber 

The  role  of  the  prechamber  is  to  precondition  the  sample  air  so 
that  transient  supersaturations  are  avoided  when  the  sample  air  enters 
the  main  chamber.  The  sample  air  is  preconditioned  approximately  to 
the  temperature  and  vapor  pressure  of  the  top  plate  (Saxena  et  al . , 
1970).  Additionally,  caution  must  be  taken  to  avoid  condensation  in 
the  sample  air  when  the  sample  enters  the  prechamber.  This  is  achieved 
by  regulating  the  filtered  air  from  the  postchamber  and  warming  it. 
Figure  7  shows  the  effect  of  this  warming. 

The  prechamber  shares  the  top  plate  with  the  other  sections.  The 
top  plate  is  constructed  of  1/8  inch  (3.2  mm)  copper  plate  with  eight 
copper  cooling  pipes  soldered  on  the  top.  The  bottom  plate  is  also 
1/8  inch  (3.2  mm)  copper  with  a  copper  cooling  coil  attached  to  the 
bottom.  The  bottom  plate  is  thermally  insulated  from  the  intake  device 
by  1/2  inch  (12.7  mm)  plexiglass  and  from  the  main  chamber  by  silicone 
sealant  of  5  mm  thickness.  The  side  walls  are  constructed  of  3.8  inch 
(9.5  mm)  plexiglass  with  the  right  side  kept  transparent  and  the  left 
side  painted  with  flat  black.  The  side  walls  are  continuous  through 
the  remainder  of  the  chamber.  The  walls  are  9  mm  high  which,  after 
2  mm  of  ice  is  applied  to  the  bottom  plate  and  1  mm  to  the  top,  give 


Figure  7.  The  effect  of  wanning  the  recycled  filtered  air. 

(a)  illustrates-  the  use  of  filtered  air  at  -10°C  while 

(b)  illustrates  the  effect  of  warming  it  to  0°C.  The 
distances  shown  are  from  the  entrance  of  the  prechamber. 
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a  6  mm  working  height  between  the  top  and  bottom  ice  layers.  The 
15  cm  length  of  the  prechamber  allows  the  sample  to  reach  steady  state 
before  entering  the  main  chamber. 

3.2.3  Main  chamber 

A  range  of  environmental  conditions  are  maintained  in  the  main 
chamber  for  activation  of  ice  nuclei.  The  shared  top  plate  of  the 
main  chamber  is  held  isothermal  while  a  temperature  gradient  is  main¬ 
tained  across  the  bottom  plate.  The  result  is  a  range  of  supersatura¬ 
tion  across  the  horizontal  median  plane  of  the  chamber  with  near 
isothermal  condition. 

The  walls  and  top  plate  of  the  main  chamber  are  the  same  as  those 
in  the  prechamber.  The  bottom  plate  is  1/2  inch  (12.7  mm)  thick  copper 
and  extends  3  cm  beyond  the  left  wall  and  4  cm  beyond  the  right  wall. 
The  use  of  a  thick  plate  at  the  bottom  enables  a  near  linear  tempera¬ 
ture  gradient  to  be  maintained  across  it.  The  bottom  plate  is  insula¬ 
ted  from  the  adjacent  plates  by  a  5  mm  width  of  silicone  sealant.  An 
angled  mirror,  placed  just  outside  the  right  wall,  allows  observation 
of  the  interior  of  the  main  chamber  while  the  chamber  is  in  operation. 

Heat  pipes  are  utilized  to  insure  temperature  uniformity  along 
the  bottom  plate  in  the  direction  of  the  flow  (Dunn  and  Reay,  1978).  In 
the  heat  pipes,  methylene  dichloride  is  used.  The  heat  pipes  are  soldered 
to  the  edges  of  the  bottom  plate.  The  top  surfaces  of  the  heat  pipes 
are  cooled  with  thermoelectric  modules  which  utilize  antifreeze  liquid 
from  the  circulating  bath  as  a  heat  sink.  The  details  of  the  heat  pipe 
and  thermoelectric  modules  installation  will  be  discussed  later. 

The  length  of  the  main  chamber  is  20  cm.  For  reasonable  flow 
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rates,  this  length  is  adequate  for  the  sample  air  to  reach  steady  state 
before  exiting  the  main  chamber.  In  most  cases,  steady  state  is  reached 
for  temperature  and  vapor  fields  in  the  first  half  of  the  chamber. 

Hence,  steady  state  conditions  are  held  for  the  last  half.  A  cross- 
sectional  drawing  is  presented  in  Figure  8. 

3.2.4  Postchamber  and  exit  device 

The  nucleated  ice  crystals  are  collected  for  counting  in  the  post¬ 
chamber.  The  new  collection  scheme  utilized  will  be  discussed  in  de¬ 
tail  later.  The  collection  scheme  however  requires  the  postchamber  to 
be  operated  at  warmer  temperatures  than  the  main  chamber  in  order  to 
avoid  spurious  ice  nucleation  on  the  ice  crystal  detection  surface, 
i.e.,  a  special  film.  Special  precautions  must  be  taken  to  avoid 
transient  supersaturations  since  the  incoming  air  from  the  main  chamber 
is  not  preconditioned  to  the  temperature  and  vapor  pressure  of  the  top 
plate.  In  order  to  avoid  these  transient  supersaturations  a  bare  metal 
surface  is  employed  on  the  bottom  plate  at  the  beginning  of  the  post¬ 
chamber.  In  this  manner,  the  heat  conduction  from  the  warmer  bottom 
plate  begins  upward  into  the  chamber  air  at  the  beginning  of  the  post¬ 
chamber  while  the  water  vapor  diffusion  process  is  delayed.  The 
height  of  the  postchamber  is  maintained  at  6  mm  for  the  first  5  cm  of 
the  postchamber  before  expanding  to  12  mm  over  the  collection  film. 

This  narrow  height  at  the  beginning  allows  the  air  to  acquire  quickly 
the  temperature  and  vapor  pressure  profiles  of  the  postchamber  and 
thereby  helps  insure  thermal  stability  of  the  air  in  the  postchamber. 

The  postchamber  height  expansion  slows  the  flow  and  allows  time  for  the 
crystals  to  settle  on  the  collecting  film. 


Figure  8.  Cross  section  of  the  main  chamber  with  the  flow  direction  ou 


Figure  9.  Schematic  of  the  flow  routing  through  the  chamber  and  support  equipment 
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The  side  walls  of  the  postchamber  are  continuous  from  the  main 
chamber.  The  bottom  plate  is  1/8  inch  (3.2  mm)  copper  with  a  cooling 
coil  on  the  bottom.  A  5  mm  high  by  5  cm  wide  copper  block  is  soldered 
to  the  bottom  plate  at  the  beginning  of  the  postchamber  to  achieve  the 
narrow  height  at  the  beginning. 

Air  is  removed  from  the  postchamber  through  an  exit  device  attach¬ 
ed  to  the  downstream  end  of  the  chamber.  It  is  constructed  of  1/8  inch 
(0.32  mm) plexiglass  and  the  design  allows  observation  into  the  chamber 
from  the  rear.  The  air  exits  through  two  pipes  in  the  bottom,  and  a 
1.5  cm  layer  of  open  cell  foam  equalizes  the  pressure  gradient  for  uni¬ 
form  flow. 

The  filtered  air  used  to  sandwich  the  sample  air  is  passed  through 
an  absolute  filter  which  consists  of  a  plexiglass  cylinder  45  cm  long 
und  2  inches  (5.1  cm)  in  diameter  with  tightly  packed  glass  wool. 

Since  the  vapor  pressure  of  the  filtered  postchamber  air  is  known,  when 
it  is  recycled  to  sandwich  the  sample  air,  the  transient  supersatura¬ 
tion  in  the  sample  air  can  be  computed  from  the  vapor  pressure  value 
and  accordingly  controlled. 

As  stated  earlier,  flow  control  and  measurement  is  crucial  for  the 
ice  nuclei  measurements  in  the  present  device.  A  schematic  of  the  flow 
conditions  is  presented  in  Figure  9.  The  flow  is  measured  by  two  flow¬ 
meters  and  accordingly  controlled.  After  the  air  exits  the  chamber, 
it  flows  through  a  meter  which  measures  and  regulates  the  total  flow, 
i.e.,  both  sample  and  filtered  air.  The  flow  then  goes  through  a 
Markson  air  pump, model  6363.  As  the  air  exits  the  pump,  it  is  divided 
into  two  flows,  one  through  the  absolute  filter  for  reuse  in  the  cham¬ 
ber  and  the  other  through  a  second  flowmeter  into  the  environment. 
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By  adjusting  the  flow  through  the  second  meter,  the  amount  of  sample 
introduced  into  the  chamber  is  regulated.  For  example,  when  the  second 
flowmeter  is  closed,  the  same  air  is  recycled  through  the  chamber  thus 
creating  a  closed  system  with  only  filtered  air  entering  the  chamber. 
When  the  second  meter  is  adjusted  to  a  flow  rate,  say  one  liter  per 
minute,  air  is  bled  from  the  system  at  a  rate  of  one  liter  per  minute 
and  is  replaced  by  sample  air  through  the  intake  device  at  the  same 
rate.  In  this  manner,  the  sample  air  enters  directly  into  the  chamber 
without  passing  through  a  meter  or  a  pump,  thus  avoiding  any  modifica¬ 
tion  of  the  sample  by  these  devices  (Langer,  1971).  The  flow  meters 
are  Manostat  and  were  calibrated  with  a  Bubble-o-Meter  in  our  labora¬ 
tory. 

3.2.5  Ice  crystal  detection 

The  detection  and  counting  of  the  nucleated  ice  crystals  in  the 
ice  thermal  diffusion  chambers  has  long  been  a  problem  (Alofs  and 
Carstens,  1976).  For  thermal  diffusion  chambers,  some  types  of  optical 
or  photographic  detection  scheme  including  the  visual  method  using 
lasers  or  other  light  sources  were  employed.  These  methods  are  awkward 
and  often  unreliable  even  for  studying  artificial  ice  nuclei  where  the 
number  of  ice  crystals  formed  can  be  controlled.  The  scarcity  of  natur¬ 
al  ice  nuclei  eliminates  these  approaches  entirely  for  their  detection. 
Hence,  a  scheme  which  collects  the  ice  crystals  is  required  for  accurate 
counting.  In  addition,  a  collection  scheme  which  provides  a  number 
integration  is  preferred,  especially  for  natural  ice  nuclei  studies, 
i.e.,  a  scheme  where  crystals  nucleated  under  the  same  environmental 
conditions  are  collected  in  the  same  region  for  counting. 
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There  are  two  critical  requirements  whicn  must  be  met  by  ary  ice 
crystal  collection  scheme,  that  all  nucleated  crystals  be  received  on 
the  collection  surface  and  that  no  additional  crystals  are  produced  by 


the  collection  surface.  It  was  found  that  an  acceptable  collection 
surface  under  the  temperature  of  the  bottom  plate  in  the  main  chamber 
was  a  difficult  task.  A  spectrum  of  materials  were  tested,  but  all  had 
significant  but  different  problems.  A  variety  of  semipermeable  mem¬ 
branes  were  tested  being  placed  on  the  ice-coated  bottom  plate  of  tne 
main  chamber  including  polyurethane,  cellophane,  cellulose  acetate  and 
cuperphane  (used  in  artificial  kidney  machines).  Most  were  acceptable 
at  warm  temperatures  but  nucleated  ice  crystals  on  the  surface  at  lower 
temperatures.  Many  fabrics  were  tested  including  silk,  velvet,  rayon, 
polyester,  satin,  cotton,  and  others.  Even  the  closest  weave  was  too 
coarse  and  allows  crystals  to  either  fall  through  or  become  obscured  by 
the  fabric.  Tissue  paper  and  carbon  paper  wrinkled  badly  wnenever  an 
ice  supersaturation  existed.  Carbon  black  dispersed  in  ether  was  poured 
over  the  ice  surface  and  after  the  ether  evaporated,  a  carbon  layer  re¬ 
mained  on  the  ice  surface.  This  method  however,  besides  being  awkward 
to  apply,  did  not  leave  a  uniform  layer  of  carbon  black. 

Mylar  copy  film  was  found  to  have  many  favorable  characteristics. 
Until  very  cold  temperatures  are  reached,  the  film  does  not  nucleate 
ice  crystals  on  the  surface.  The  black  surface  offers  a  high  contrast 
with  the  ice  crystals  under  low  angle  illumination.  The  backing  is 
nonpermeable  with  water  vapor,  a  characteristic  used  to  our  advantage 
in  the  collection  scheme.  The  film  is  hardy,  readily  available,  and 
requires  no  special  handling  or  preparation.  The  surface  is  smooth  and 
does  not  wrinkle  as  does  normal  carbon  paper.  In  addition,  it  was 
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found  that  the  surface  has  a  poisoning  effect  on  the  nuclei  after  the 
collected  crystals  melted  and  evaporated.  This  allows  the  film  to  be 
reused  many  times  without  fear  of  background  crystal  counts. 

The  film  is  adhered  to  a  1/32  inch  (0.8  mm)  aluminum  sheet  which 
can  be  easily  slid  into  the  postchamber  from  the  rear.  The  film  is  re¬ 
moved  from  the  postchamber  after  the  crystals  have  been  collected  and 
put  in  a  freezer  for  photographing  of  the  collected  ice  crystals. 

Figure  10  shows  the  method  of  photographing.  The  film  is  then  allowed 
to  dry  for  reuse. 

3.2.6  Temperature  control  and  measurement 

A  Neslab  circulating  refrigerated  bath.  Model  RTE-3,  furnishes  the 
temperature  control  for  the  top  plate  and  the  bottom  plates  of  the  pre¬ 
chamber  and  postchamber.  Furthermore,  it  is  used  as  a  heat  sink  for 
the  thermoelectric  modules.  The  coolant,  commercial  antifreeze,  is 
routed  from  the  circulating  bath,  first  through  the  coils  on  the  bottom 
plates,  then  across  the  thermoelectric  modules  where  heat  is  gained, 
and  finally  through  the  coolant  pipes  on  the  top  plate.  This  routing 
scheme  permits  the  top  plate  to  stay  at  a  temperature  approximately  one 
degree  warmer  than  the  bottom  plates  of  the  prechamber  and  postchamber 
and  thereby  enhances  the  thermal  stability  of  the  flow  and  secures  a 
positive  vertical  vapor  pressure  gradient  in  the  postcnamber.  In  this 
manner  sublimation  of  the  collected  ice  crystals  is  avoided  in  the  post¬ 
chamber. 

Cambion  thermoelectric  modules,  model  2001,  are  used  with  the  heat 
pipes  for  temperature  control  on  the  edges  of  the  main  chamber  bottom 
plate.  They  are  installed  using  thermal  compound  with  the  cooling 
sides  in  contact  with  the  tops  of  the  heat  pipes  and  the  warm  side  in 


Q. 

2 

< 


Figure  10.  The  arrangement  for  photographing  the  crystals  on  the  collection  film. 
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contact  with  a  plate  cooled  by  the  coolant  from  the  circulating  bath. 
Four  modules  are  used  for  the  cold  side  and  two  for  the  warm  side.  The 
modules  are  wired  in  series  and  their  output  is  controlled  by  adjusting 
the  current  through  them.  Separate  DC  power  supplies  are  used  for  each 
side.  In  addition  to  the  cooling  modules,  a  heating  wire  is  attached 
to  the  edge  of  the  warm  side  to  supply  heat  when  required.  The  wire  is 
a  40  cm  length  of  24  gauge  constantan.  The  same  power  supply  used  for 
the  thermoelectric  modules  on  that  side  is  used  to  heat  the  wire. 

The  temperature  measurements  are  made  by  copper-constantan  thermo¬ 
couples  with  24  gauge  extention  wires.  The  thermocouples  were  spot- 
welded  to  obtain  a  single  contact  between  the  two  metals.  A  melting 
ice  bath  is  used  for  the  temperature  reference.  The  voltage  output  of 
the  thermocoupl es  is  measured  using  a  Keithley  digital  voltmeter.  Model 
174.  This  combination  produces  an  accuracy  of  0.0 1°C  at  -10°C. 

The  temperature  is  measured  at  nine  locations  in  the  chamber.  The 
temperature  recording  is  automated  and  controlled  by  a  Terak  minicompu¬ 
ter.  Each  temperature  is  recorded  every  25  seconds  on  a  magnetic  disk 
and  shown  on  a  video  display  screen.  The  temperature  is  measured  at 
the  following  locations: 

-  Prechamber;  middle  of  the  top  and  bottom  plates. 

-  Main  chamber;  half  way  down  the  length  of  the  bottom  plate  at 

the  left  wall,  middle  of  the  chamber  and  right 
wall  and  at  the  middle  of  the  top  plate. 

-  Postchamber;  middle  of  the  bottom  plate. 

-  Filtered  air;  just  prior  to  the  entrance  into  the  intake  device. 

-  Sample  air;  just  prior  to  the  entrance  into  the  intake  device. 
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The  thermocouples  which  measure  the  temperature  in  the  plates  are 
painted  with  a  solution  of  plexiglass  to  electrically  insulate  them 
from  the  copper  plates.  In  addition,  they  were  coated  with  thermal 
compound  before  being  inserted  into  holes  in  the  plates.  The  thermo¬ 
couples  which  measure  air  temperatures  are  inserted  directly  into  the 
air  flow. 

3.3  Chamber  Operation 

3.3.1  Application  of  ice  surfaces 

The  ice  surfaces  on  the  top  and  bottom  plates  must  be  of  uniform 
thickness  and  have  smooth  surfaces.  Many  methods  for  applying  the  ice 
were  tested  but  most  of  them  produced  either  non-uniform  ice  layers  or 
ice  layer  with  rough  surfaces.  Two  different  methods  were  finally 
chosen  to  apply  the  bottom  and  top  ice  layers.  For  the  bottom  plates, 

3 

110  cm  of  water  are  poured  over  the  bottom  plates  where  ice  coating  is 
required.  This  gives  a  2  mm  layer  of  ice  over  these  plates.  Capillary 
action  of  the  water  keeps  it  from  running  off  the  front  and  near  edges. 
The  copper  plate  used  at  the  beginning  of  the  postchamber  is  placed  in 
the  water  and  frozen  in  place.  As  the  plates  are  cooled,  the  water 
freezes  in  a  uniform  layer  with  a  smooth  surface. 

For  application  of  the  ice  layer  on  the  top  plate,  the  top  plate 
is  inverted  and  tape  placed  along  the  outside  5  mm  of  each  side.  Next 
the  plate  is  cooled  to  -5°C.  Water  is  applied  to  the  surface  using 
a  sponge  covered  with  a  plastic  mesh  like  those  used  in  service  sta¬ 
tions  to  clean  automobile  windshields.  The  thin  water  layer  applied 
soon  freezes.  The  procedure  is  repeated  until  a  1  mm  layer  of  ice 
builds  up.  Then,  the  tape  is  removed  while  moving  the  wet  spcnue  along 
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the  tape  and  warming  the  ice.  This  procedure  prevents  any  chipping  of 
the  ice.  This  ice-free  area  on  the  edges  where  the  tape  is  removed, 
is  used  to  create  a  leakproof  seal  with  the  top  of  the  side  walls. 
Figures  11(a)  and  11(b)  show  this  ice  application  scheme. 

3.3.2  Sandwiching  of  the  sample 

The  main  chamber  is  supersaturated  with  respect  to  ice  throughout 
the  interior  except  on  the  boundaries  where  ice  saturation  is  main¬ 
tained.  The  ice  supersaturation,  however,  is  not  constant  vertically 
between  the  top  and  bottom  plates  but  ranges  from  zero  at  the  boundaries 
to  a  maximum  near  the  center.  For  sufficiently  large  temperature  dif¬ 
ferences  between  the  plates,  supersaturation  with  respect  to  water 
occurs  with  a  maximum  also  near  the  center  but  not  necessarily  colocated 
with  the  maximum  for  ice  supersaturation.  However,  for  a  thin  sample 
layer  around  the  center,  both  ice  and  water  supersaturation  are  nearly 
constant  at  values  close  to  their  respective  maxima. 

The  temperature  varies  almost  linearly  between  the  plates  of  the 
main  chamber.  Hence,  within  the  thin  sample  layer,  the  temperature 
varies  little.  The  thinner  the  layer,  the  less  the  variance  of  the 
temperature  through  the  layer.  Therefore  the  sample  should  be  confined 
to  as  shallow  a  layer  as  possible  for  near  isothermal  conditions.  How¬ 
ever,  there  is  a  practical  limit  on  the  thickness  of  the  layer  since, 
as  the  sample  layer  thickness  is  reduced,  the  volume  sampling  rate  be¬ 
comes  reduced.  Thus,  for  a  given  sampling  rate,  there  is  a  correspond¬ 
ing  minimum  sample  layer  thickness. 

Under  the  shear  flow  condition  of  the  chamber,  the  velocity  varies 


from  zero  at  the  boundaries  to  a  maximum  at  the  vertical  center.  At 
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Figure  11(a).  Application  of  the  ice  layer  on  the  top  plate  at  -5°C. 

Figure  11(b).  Removal  of  the  tape  along  the  edge  of  the  top  plate 

where  an  ice-free  surface  is  required  for  contact  with 
the  side  wall . 
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the  center  of  the  flow  there  exist  a  thin  region  where  the  flow  veloc¬ 
ity  is  nearly  constant  at  the  maximum  velocity  rate.  Therefore,  if  the 
sample  is  confined  to  this  thin  region  near  the  center  of  the  chamber, 
it  will  be  exposed  to  nearly  uniform  fields  of  supersaturation,  temper¬ 
ature,  and  velocity.  As  discussed  earlier,  the  temperature  and  water 
vapor  fields  and  therefore  the  supersaturation  can  be  accurately  con¬ 
trolled  and  the  velocity  field  can  be  computed. 

The  correct  positioning  and  depth  of  the  sample  layer  is  done  by 
means  of  the  intake  device  described  earlier.  By  controlling  the  rela¬ 
tive  thickness  of  the  filtered  air  layers,  the  vertical  position  of  the 
sample  layer  between  the  plates  can  be  adjusted.  The  thickness  of  the 
sample  layer  is  controlled  by  adjusting  the  bleed  air  through  the  sec¬ 
ond  flow  meter  so  that  the  desired  sampling  rate  is  achieved.  The 
vertical  position  of  the  sample  layer  is  adjusted  to  the  region  where 
the  supersaturation,  temperature,  and  velocity  fields  are  nearly  con¬ 
stant.  Figure  12  illustrates  the  profiles  of  these  parameters  and  the 
positioning  of  the  sample.  The  intake  device  also  introduces  the  fil¬ 
tered  air  between  the  sample  air  and  the  wall  so  that  wall  effects  are 
avoided. 

3.3.3  Preconditioning  of  the  sample 

Caution  must  be  taken  to  insure  that  condensation  does  not  occur 
in  the  sample  air  prior  to  entering  the  main  chamber  since  condensation 
could  alter  the  number  of  ice  nuclei  in  the  sample  or  modify  their  be¬ 
havior  by  providing  an  uncontrolled  opportunity  for  the  liquid  phase 
to  build  up  on  them.  Since  the  temperatures  of  the  filtered  and  the 
sample  airs  are  measured,  if  the  vapor  pressure  of  the  sample  were  known, 
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computations  using  the  diffusion  equation  with  advection  could  determine 
if  water  saturation  were  ever  reached.  Instead  of  using  an  additional 
device  for  measuring  the  vapor  pressure  of  the  sample  air,  we  assume  a 
water  saturated  sample  and  proceed  with  the  calculations.  Under  condi¬ 
tions  where  transient  water  supersaturations  are  predicted  to  occur,  the 
filtered  air  is  heated  properly  to  avoid  them.  By  heating  the  filtered 
air,  the  temperature  field  of  the  predried  filtered  air  remains  suffi¬ 
ciently  high  while  the  vapor  diffuses  into  it  so  that  water  supersatura¬ 
tion  does  not  occur.  This  effect  was  previously  shown  in  Figure  7. 

The  purpose  of  the  prechamber  is  to  precondition  the  sample  to  a 
temperature  and  a  vapor  pressure  close  to  those  of  the  top  plate.  This 
preconditioning  prevents  transient  supersaturations  from  happening  in 
the  sample  air  when  entering  into  the  main  chamber.  The  final  condition 
attained  in  the  prechamber  is  slightly  different  from  isothermal  since 
the  top  plate  is  held  about  1°C  warmer  than  the  bottom.  However,  our 
computations  show  that  this  preconditioned  profile  is  adequate  to  avoid 
transient  supersaturations. 

3.3.4  Activation  of  ice  nuclei 

The  sample  moves  through  a  region  in  the  main  chamber  where  the 
supersaturation  and  temperature  are  known.  The  residence  time  is  suffi¬ 
cient  for  the  ice  supersaturation  value  of  the  sample  air  to  come  within 
10%  of  the  steady  state  value.  In  most  cases,  this  is  accomplished  in 
the  first  one  half  of  the  main  chamber  length.  Therefore,  the  sample 
air  stays  within  10%  from  the  steady  state  value  of  ice  supersaturation 
for  the  last  half  of  the  main  chamber.  It  is  in  this  region  where  ice 
nucleation  occurs  in  the  sample  air.  With  the  horizontal  variation  of 


supersaturation  across  the  chamber,  the  ice  nucleations  which  occur 
at  different  positions  across  the  chamber  have  correspondingly  differ¬ 
ent  supersaturation  thresholds. 

The  nucleated  crystals  are  carried  out  of  the  main  chamber  by  the 
flow  and  settle  on  the  collection  film  in  the  postchamber.  We  stratify 
the  horizontal  supersaturation  variation  by  defining  lanes  of  nearly 
constant  supersaturations  horizontally  across  the  flow.  The  lanes  ex¬ 
tend  down  the  main  chamber  and  into  the  postchamber.  Because  laminar 
flow  is  maintained,  the  crystals  collected  on  the  film  in  each  lane 
must  have  nucleated  at  the  supersaturation  and  temperature  of  that  lane 
in  the  main  chamber. 

3.3.5  Collection  of  the  ice  crystals  and  analysis 

The  purpose  of  the  postchamber  is  to  collect  the  ice  crystals 
which  nucleated  in  the  main  chamber.  There  are  two  levels  of  the  bot¬ 
tom  plates  in  the  post  chamber-,  the  first  maintains  the  same  level  be¬ 
tween  ice  layers  as  the  main  chamber  and  the  second  is  lower  so  that 
the  flow  slows  down  there  and  crystals  have  time  to  settle  on  the  film 
before  the  sample  air  exits  the  chamber.  The  temperature  profile  in 
the  postchamber  is  the  same  as  that  in  the  prechamber,  i.e.  the  bottom 
plate  is  maintained  about  1°C  colder  than  the  top  plate.  This  profile 
increases  the  thermal  stability  of  the  air  and  maintains  a  downward 
flux  of  vapor  so  that  the  collected  crystals  do  not  sublimate.  This 
arrangement  may  lead  to  problems;  as  the  cold  air  from  the  main  chamber 
moves  over  the  warm  bottom  plate  of  the  postchamber,  thermally  unstable 
conditions  can  result.  Furthermore,  the  air  adjacent  to  the  bottom 
plate  in  the  main  chamber  has  a  very  low  vapor  pressure  and  when  it 
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moves  to  the  postchamber,  it  can  sublimate  the  collected  crystals  un¬ 
less  the  vapor  pressure  is  increased.  This  cold  dry  air  can  also  eva¬ 
porate  or  freeze  the  water  droplets  deposited  on  the  collection  film. 
For  these  reasons,  the  first  5  cm  of  postchamber  are  maintained  at  the 
6  mm  height  with  both  the  top  and  bottom  plates  covered  with  ice.  This 
small  height  suppresses  the  thermal  instability  by  making  the  tempera¬ 
ture  field  reach  the  new  steady  state  more  quickly  and  reducing  the 
available  time  for  thermal  instability.  It  also  permits  the 
vapor  field  of  the  air  to  reach  the  postchamber  profile.  The  length 
of  this  ice  covered  bottom  surface  must  sufficiently  be  short,  however, 
so  that  nucleated  crystals  fall  on  the  collection  film  instead  on  the 
ice  surface. 

The  movement  of  this  supersaturated  cold  air  into  the  relatively 
warm  postchamber  will  induce  undesirable  transient  supersations  unless 
precautions  are  taken  (Fitzgerald,  1970;  Saxena  et  al ,  1970).  These 
transient  conditions  occur  because  the  vapor  diffuses  faster  than  the 
heat  conducts.  To  avoid  this  problem,  the  vapor  diffusion  is  deliber¬ 
ately  delayed  at  the  beginning  of  the  postchamber  while  heat  conduction 
proceeds  in  the  air.  This  delay  is  acheived  by  substituting  bare  (non¬ 
ice  covered)  a  copper  strip  for  the  ice  layer  on  the  bottom  plate  at 
the  beginning  of  the  postchamber.  Computations  show  this  procedure  is 
adequate  for  avoiding  supersaturations. 

The  Mylar  copy  film  for  ice  crystal  collection  is  used  by  adhering 
it  to  a  1/32  in  aluminum  plate  with  Scotch  double-coated  tape.  It  is 
introduced  into  the  postchamber  by  sliding  from  the  rear.  When  the 
temperature  of  the  postchamber  is  above  -10°C,  best  results  are  obtain¬ 
ed  by  using  a  pretreated  Mylar  copy  film  which  carries  condensed  water 
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droplets.  The  pretreatment  procedure  is  as  follows;  the  Mylar  film  is 
placed  under  a  plexiglass  cover  which  has  wetted  filter  paper  adhered 
to  the  interior.  The  cover  is  kept  at  room  temperature.  The  film  and 
the  cover  are  placed  on  a  thick  metal  plate  precooled  in  a  freezer. 

The  film  cools  first  and  vapor  diffuses  to  the  film  surface  from  the 
filter  paper  of  the  cover.  The  film  becomes  coated  with  small  super¬ 
cooled  water  droplets.  When  the  postchamber  temperature  is  below  -10°C, 
the  film  should  be  used  without  pretreatment. 

In  order  to  confirm  that  the  Mylar  film  produces  no  crystals  on 
its  surface,  blank  tests  were  run  using  only  filtered  air  in  the  cham¬ 
ber.  The  tests  consistently  produced  no  background  ice  crystals  on 
the  film.  This  confirmed  that  the  ice  crystals  collected  during  ice 
nuclei  sample  tests  were  indeed  a  result  of  ice  nucleation  in  the  main 
chamber. 

An  Olympus  35  ran  0M-1  camera  was  used  for  recording  ice  crystals 
on  the  Mylar  film  and  their  subsequent  analysis.  After  a  sufficient 
sample  volume  is  drawn  through  the  chamber  and  the  ice  nuclei  are  acti¬ 
vated,  the  film  is  placed  under  a  thin  plexiglass  cover  whose  interior 
is  coated  with  ice  and  is  then  moved  to  a  freezer  for  photographing. 

The  plexiglass  cover  protects  ice  crystals  on  the  film  while  transferring 
to  the  freezer.  The  film  is  placed  on  a  metal  plate  in  the  freezer  held 
at  about  -10°C.  If  the  crystals  are  not  large  enough  to  photograph, 
the  film  is  placed  under  a  plexiglass  cover  1  cm  high  with  wetted  fil¬ 
ter  paper  adhered  to  the  interior  to  grow  the  crystals  to  easily  recog¬ 
nizable  sizes.  The  cover  is  kept  at  room  temperature.  When  the  crys¬ 
tals  become  sufficiently  large,  they  are  illuminated  with  a  lamp  at  a 
low  angle.  The  crystals  are  photographed  at  about  a  70  degree  angle 
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from  the  plane  of  the  film.  The  photographing  arrangement  was  previous¬ 
ly  shown  in  Figure  10. 

After  the  photographs  are  developed,  the  lanes  and  their  supersa¬ 
turations  are  identified  using  an  overlay.  The  number  of  crystals  in 
each  lane  is  determined  and  recorded.  These  number  counts  combine  with 
the  computed  temperature  and  supersaturation  fields  are  used  for  data 
analysis.  An  example  of  the  collected  crystals  is  presented  in  Figure 
13(a)  and  the  use  of  the  overlay  shown  in  Figure  13(b). 

Another  source  of  disturbance  in  the  flow  was  dendrites  grown  on 
the  bottom  plate  of  the  main  chamber.  This  dendrite  growth  is  not  un¬ 
expected  since  a  continuous  downward  flux  of  vapor  exist  as  well  as  ice 
supersaturation  above  the  bottom  plate.  These  crystals  protrude  into 
the  chamber  interior  and  can  disturb  the  laminar  flow.  These  distur¬ 
bances  were  confirmed  with  the  smoke  tests. 

Two  precautions  are  taken  to  eliminate  this  dendrite  growth. 

First,  whenever  the  chamber  is  cooled  to  operating  temperatures  and  not 
in  use,  a  plastic  sheet  is  inserted  between  the  plates.  The  dendrites 
then  grow  on  the  sheet  and  are  removed  from  the  chamber  with  the  sheet 
before  measurements.  However,  while  the  chamber  is  used  for  sampling, 
crystals  do  begin  to  grow  on  the  bottom  plate.  These  crystals  are 
melted  periodically  by  moving  a  warm  copper  plate  over  the  surface. 

The  plate  is  thermally  insulated  from  the  top  ice  surface  completely 
and  only  exposed  to  the  bottom  surface  for  1  cm  of  its  10  cm  length. 

This  design  permits  good  control  of  the  dendrite  melting  process.  The 
short  exposed  length  of  the  copper  plate  permits  a  quick  melting  and  re¬ 
freezing  of  the  crystals,  thereby  leaving  a  smooth  surface  after  melt¬ 
ing  no  more  ice  layer  surface  than  necessary. 
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Figure  13(a).  Example  of  the  received  ice  crystals  on  the  collection 
film.  The  sample  was  silver  iodide. 

Figure  13(b).  Use  of  the  overlay  dividing  the  collection  film  into 
lanes.  The  number  of  collected  ice  crystals  in  each 
lane  is  counted. 
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The  shear  flow  velocity  profile  resulting  from  the  presence  of 
the  top  and  bottom  plates  has  been  previously  discussed.  Shear  flow 
also  exists  along  the  side  walls,  although  its  depth  of  penetration  into 
the  chamber  is  small.  This  velocity  shear  induces  significant  thermal 
instability  along  the  walls.  The  air  adjacent  to  the  wall  moves  more 
slowly  than  that  of  the  interior  and  therefore  requires  a  shorter  dis¬ 
tance  down  the  chamber  to  stabilize.  Consequently,  in  the  prechamber 
and  main  chamber,  the  air  close  to  the  wall  becomes  colder  than  that 
of  the  inner  flows  and  causes  subsidence  along  the  wall.  In  the  post¬ 
chamber,  the  situation  is  vice  versa;  air  close  to  the  wall  warms 
faster  and  convection  occurs  at  the  wall.  The  height  of  the  main 
chamber  was  originally  built  to  be  1  cm  and  that  of  the  postchamber 
2  cm.  For  these  heights,  the  thermal  instabilities  along  the  walls 
were  very  large.  In  the  postchamber ,  the  air  adjacent  to  the  wall  at 
the  entrance  was  carried  halfway  across  the  top  plate  by  the  time  it 
exited  the  chamber.  A  photograph  of  a  smoke  test  showing  this  insta¬ 
bility  is  given  in  Figure  14(a). 

The  solution  to  this  problem  was  twofold.  First,  the  heights  were 
lowered  to  6  mm  and  12  mm,  respectively,  in  the  main  chamber  and  post- 
chamber  so  as  to  reduce  the  subsidence  or  convective  force  as  well  as 
its  duration.  This  alone  eliminated  most  of  the  instability  but 
as  can  be  seen  in  Figures  14(b)  and  14(c),  sample  air  at  the  wall 
is  subject  to  some  remaining  convective  forces.  However,  the 
width  of  the  filtered  air  separating  the  sample  from  the  walls  was 
also  increased  slightly  to  insure  that  the  sample  air  stayed  out  of 
any  additional  instability.  These  precautions  proved  adequate  for  main¬ 
taining  laminar  flow  in  the  sample  layer  as  was  previously  shown  in 
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Figure  14(a).  Photograph  of  a  smoke  test  where  the  sample  layer 

extended  to  the  wall.  A  plexiglass  plate  is  substituted 
for  the  top  plate  and  the  chamber  is  viewed  from  above. 
The  air  flow  is  from  left  to  right. 

Figure  14(b)  and  (c).  Photograph  of  a  smoke  test  viewed  from  the 
exit  end  of  the  chamber.  The  sample  layer  extended 
to  the  wall.  The  smoke  is  illuminated  at  the  entrance 
to  the  postchamber  in  (b)  and  at  the  end  of  the  post¬ 
chamber  in  (c). 
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Figure  9. 

A  more  obvious  source  of  turbulence  would  be  leaks  in  the  chamber. 
All  parts  which  are  not  routinely  removed  were  sealed  when  they  were 
installed  with  either  foam  strips  or  silicone  rubber  sealant.  The 
top  plate  and  the  exit  device  are  routinely  removed.  The  top  plate  is 
sealed  by  applying  a  water  bead  along  the  edges.  This  bead  quickly 
freezes  to  form  an  ice  seal.  Closed  cell  foam  weather  stripping  is 
used  to  seal  between  the  exit  device  and  the  postchamber.  Smoke  tests 
confirm  the  leakproof  sealing  of  the  chamber  and  the  laminar  flow 
throughout  the  chamber  as  shown  in  Figure  15(a)  and  15(b). 

3.4  Wall  Effect  Determination 

Effects  of  side  walls  on  the  temperature  and  vapor  fields,  and 
therefore  the  supersaturation  fields  have  long  been  recognized  as  a 
significant  problem  in  thermal  diffusion  chambers  (Elliott,  1971),  but 
only  recently  have  quantitative  calculations  been  carried  out  to  deter¬ 
mine  the  penetration  depth  of  the  effects  into  the  interior  of  the 
chambers.  These  wall  effects  are  attributed  to  the  ice  saturated  con¬ 
dition  on  the  walls  due  to  formed  ice  crystals  on  them,  i.e.  the  near 
linear  vapor  pressure  profile  of  the  interior  curves  to  the  saturation 
values  on  the  walls.  Computations  similiar  to  those  presented  in  the 
paper  by  Tomlinson  and  Fukuta  (1979)  were  carried  out  for  the  rectangu¬ 
lar  geometry  of  the  present  ice  thermal  diffusion  chamber.  The  wall 
effects  were  included  in  every  computation  of  the  temperature  and 
supersaturation  fields,  i.e.  for  each  sample  taken,  the  complete  calcu¬ 
lation  was  made  using  the  measured  temperatures.  In  this  manner, 
accurate  values  for  the  temperature  and  supersaturation  in  the  region 
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Figure  15.  Photograph  of  a  smoke  test  where  the  top  plate  was 

replaced  with  plexiglass.  In  (a),  pressure  waves  were 
created  by  tapping  the  top  of  the  intake  device  and  the 
wave  fronts  were  observed  as  they  were  advected  down¬ 
stream  by  the  flow.  In  (b),  smoke  was  released  through 
pinholes  in  a  straw  across  the  entrance  to  the  prechamber 
and  the  flow  pattern  of  the  smoke  lines  was  observed. 
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where  ice  nucleation  occurs  are  obtained.  Appendix  A  includes  the 
equations  and  method  used  for  these  computations. 

3.5  Error  Analysis 

The  errors  which  are  present  in  the  chamber  are  due  primarily  to 
those  in  temperature  measurements.  As  mentioned  earlier,  the  error 
in  temperature  measurement  with  the  equipment  used  is  0.01°C  at  -10°C. 
The  temperature  of  the  chamber  was  normally  very  stable  when  cooled 
to  steady  state.  Small  drifts  were  sometimes  observed  to  occur,  but 
were  limited  to  ±0.02°C.  For  the  error  analysis,  the  value  of  ±0.03°C 
is  used.  Trial  runs  were  made  using  the  extremes  of  the  temperature 
errors.  The  following  ranges  of  relative  error  were  determined: 

18  -  30  %  for  0.0  ^  S  {%)  ^  1 .0 

w 

10  -  18  %  for  1 .0  -  $.{%)  ^  5.0 

W 

1  -  2  %  for  12.0  *  s  (%)  *  25.0 

The  errors  associated  with  the  wall  effect  have  been  determined 
and  accordingly  the  sample  air  removed  from  the  affected  region.  There¬ 
fore,  there  is  no  wall  effect  error  in  the  sample.  The  computational 
errors  are  small  compared  with  the  above  ranges  of  error  and  are 
negl igible. 


CHAPTER  4 


NUMERICAL  MODELING  OF  THE  CHAMBER  OPERATION 

A  detailed  knowledge  of  the  conditions  throughout  the  chamber  in¬ 
terior  is  essential  in  order  for  any  meaningful  detailed  studies  to  be 
carried  out.  Numerical  models  which  simulate  the  conditions  can  furnish 
the  needed  knowledge  if  the  boundary  conditions  are  properly  described 
under  steady  state.  Two  such  models  have  been  developed  for  use  in 
this  study.  The  first  computes  the  steady  state  conditions  in  vertical 
planes  perpendicular  to  the  flow  direction.  The  other  computes  the 
steady  state  conditions  in  vertical  planes  parallel  to  the  flow  direc¬ 
tion.  The  mathematical  approaches  used  in  these  models  are  quite  dif¬ 
ferent  but  their  computed  fields  agree  very  well  along  lines  where  the 
planes  of  each  model  intersect,  indicating  consistency  of  both  models. 

4.1  Model  Describing  the  Chamber  Conditions  under 

Steady  State  and  Wall  Effects 

The  first  numerical  model  was  designed  to  compute  the  steady  state 
conditions  in  a  plane  perpendicular  to  the  flow  with  considerations 
given  to  wall  effects.  Twomey  (1963)  concluded  that  the  influence  of 
the  side  walls  on  the  interior  fields  of  the  chamber  must  be  considered 
when  computing  the  fields.  The  present  approach  is  similiar  to  that 
presented  by  Tomlinson  and  Fukuta  (1979)  but  employs  rectangular  geome¬ 
try  for  the  chamber.  As  discussed  in  the  latter  paper,  the  models 
assumes  wettable  walls  made  of  a  thermal  conductor.  Vapor  pressures 
values  were  used  instead  of  vapor  density  for  better  accuracy  as 
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suggested  by  Katz  and  Mirabel  (1975). 

For  the  steady-state  transportation  of  water  vapor  without  sources 

2 

or  sinks,  the  Laplace  equation  for  vapor  pressure  7  P  =  0  holds  with 
sufficient  accuracy  (Fukuta  and  Saxena,  1979).  The  approach  for  T  and 
P  computations  is  as  follows:  Since  P  is  of  an  additive  nature,  the 
solution  of  the  Laplace  equation  for  P  under  the  original  boundary 
conditions  is  broken  down  into  five  component  pressures  satisfying 
boundary  conditions  (1)  through  (5)  in  Figure  16  (Carslaw  and  Jaeger, 
1959)  and  the  Laplace  equation  for  each  of  them  is  solved  separately 
and  the  solutions  are  added  later.  On  the  left  wall. 


PL(z)  =  PBL  exp[^  (T^  -  jjjy)] 


where  L  is  the  latent  heat  of  sublimation  per  gram  water  and  Ry  the 
specific  gas  constant  for  water  vapor,  Tg^  the  temperature  on  the  left 
side  of  the  bottom  plate  and  TL(z)  the  vertical  temperature  profile  on 
the  left  wall.  A  similiar  equation  holds  for  the  vapor  pressure  pro¬ 
file  on  the  right  wall.  The  vertical  temperature  fields  on  both  walls 
and  therefore  interior  of  the  chamber  are  assumed  to  be  linear.  The 
equations  and  procedure  used  for  solving  for  each  component  pressure 
P^  through  P^  is  presented  in  Appendix  A. 

The  solution  for  the  Laplace  equation  for  temperature  under  a 
linearly  changing  vertical  wall  temperature  profile  is  simply  a  linear 
function  of  z.  The  saturation  vapor  pressure  Pg  is  computed  using  the 
sixth  degrees  polynomials  from  Rasmussen  (1978).  The  vapor  pressure  P 
at  any  point  inside  the  chamber  is  given  by 


P(x,z)  =  P^  +  P2(x,z)  +  P3(x,z)  +  P4(x,z)  +  ?5(x,z), 


63 


64 


Using  this  pressure  and  the  computed  saturation  vapor  pressure,  the 
supersaturation  fields  are  computed  for  the  chamber  interior.  In  turn, 
these  fields  are  used  to  determine  the  extent  of  the  wall  effects. 

This  model  was  run  for  the  ranges  of  temperatures  used  in  this 
study.  It  was  found  that  the  effects  of  the  walls  were  confined  to 
within  1.5  cm  of  the  walls,  i.e.  the  supersaturation  values  were  within 
10%  of  the  infinte  parallel  plates  value  for  all  points  further  than 
1.5  cm  from  the  walls.  Based  on  this  model,  it  was  decided  that  no 
sample  air  should  be  placed  closer  to  the  walls  than  1.5  cm. 

The  measurements  taken  from  the  chamber  during  experiments  indica¬ 
ted  that  although  a  thick  copper  plate  was  used  for  the  bottom,  devia¬ 
tions  from  a  linear  temperature  profile  were  evident  on  it.  Since  the 
temperature  is  measured  at  three  locations  across  the  plate,  a  quadratic 
function  is  used  to  describe  the  temperature  profile.  This  should  give 
a  good  approximation  to  the  actual  temperature  profile  and  the  accuracy 
of  the  model  is  preserved. 

The  model  is  run  for  every  experiment,  thereby  eliminating  the  re¬ 
quirement  that  the  temperature  be  regulated  exactly  to  some  preset  set 
of  values  for  which  the  model  has  been  previously  run.  The  temperatures 
measured  for  each  experiment  yield  computed  sets  of  fields.  The  model 
requires  approximately  25  minutes  to  run  on  the  Terak  for  each  set  of 
experimental  data  and  computes  the  following  fields: 

1)  Temperature 

2)  Vapor  pressure 

3)  Ice  supersaturation 

4)  Water  supersaturation 

5)  Error  (or  deviation)  from  the  infinite  plate  case  for  ice 


supersaturation 

6)  Error  (or  deviation)  from  the  inf inite  plate  case  for  water 
supersaturation 

Normally  only  the  temperature,  ice  supersaturation  and  water  supersa¬ 
turation  fields  are  saved  for  output.  An  example  of  the  computed  tem¬ 
perature  and  water  supersaturation  fields  is  shown  in  Figure  17. 


4.2  Model  Describing  Shear  Flow  and 

Transient  Supersaturation 

The  second  numerical  model  developed  computes  the  steady  state 
conditions  in  a  vertical  plane  parallel  to  the  flow  direction  with  con¬ 
siderations  given  to  the  shear  flow  profile,  the  temperature  and  total 
pressure  dependencies  of  the  vapor  and  thermal  diffusivities.  Previous 
computations  by  other  investigators  have  dealt  with  slab  flow  with  con¬ 
stant  values  for  the  coefficients  (Saxena  et  al.,  1970;  Fitzgerald, 
1970;  Brown  et  al . ,  1979;  Fukuta  and  Saxena,  1979a). 

The  shear  flow  velocity  profile  between  parallel  plates  is  given 


V(z) 


1  dPa 

^-d7z(H-z) 


0) 


dP 

where  y  is  the  viscosity  of  the  fluid,  the  total  pressure  gradient 
in  the  direction  of  the  flow  and  H  the  distance- between  the  plates 
(Clark  and  Kays,  1953;  Schlichting,  1955;  Bird,  1960:  White,  1974). 

The  z  direction  is  taken  in  the  vertical  between  the  plates  and  y  in  the 
direction  of  the  flow.  This  profile  is  referred  to  as  Poiseuille  flow 
and  has  a  non-slip  condition  at  the  boundaries  and  a  maximum  velocity 


midway  between  the  plates. 

The  equations  used  for  heat  and  vapor  transport,  respectively,  are 
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AT  2  AT 

~  =  o  T  +  7<VT  -  V  r1 

3 1  ' - v - J  ay 

a 

~  =.D72P  +  7D7P  -  V  —  , 

d  t  v - 1 — V — - - '  ay 

b 


(2) 

(3) 


where  k  and  0  are  the  thermal  and  vapor  diffusivities ,  respectively 
(Jakob  and  Hawkins,  1  957  ;  Kays,  1966;  Kutateladze,  1966;  Gebhart,  1971  j. 
Scale  analysis  indicate  that  the  horizontal  components  of  the  tern.;  in 
(a)  of  Eq.  (2)  and  (b)  of  Eq.  (3)  are  negligible  when  compared  to  the 
vertical  components.  Therefore,  Eqs.  (2)  and  (3)  reduce  to 


3T  _  afj  3k 

3t  <  -2'  3z 

o  Z 

3P  .  n  32P  30 

3t  '  U  .,2  3z 

aZ 


(4) 

(5) 


The  resulting  equations  are  second  order  non-linear  partial  differ¬ 
ential  equations.  Since  what  is  needed  is  the  steady  state  profiles  for 
the  interior  of  the  chamber  excluding  the  initial  time  dependent  parts, 
setting  3 T / 5 1  =  0  in  Eq.  (4)  and  3 P/3 1  =  0  in  Eq.  (5)  gives  the  steady 
state  profiles.  Equations  (4)  and  (5)  were  integrated  numerically  under 
the  measured  boundary  conditions  of  the  experiments.  The  numerical 
schemes  used  are  presented  in  Appendix  B.  The  integration  schemes  are 
stable  and  converge  nicely.  The  minimum  simulated  time  used  was  that 
required  for  the  slowest  parcel  to  be  advected  through  the  chamber. 

This  turned  out  to  be  always  longer  than  the  time  required  to  reach 
steady  state  if  the  chamber  had  been  static. 

The  primary  uses  for  this  model  are  to  determine  exactly  what  the 


68 


transient  conditions  are  and  where  steady  state  is  reached  in  each 
section  of  the  chamber.  Furthermore,  the  model  was  used  to  determine 
to  what  temperature  the  recycled  filtered  air  needed  to  be  warmed  in 
order  to  avoid  condensation  in  the  prechamber  as  shown  in  Figure  7  of 
Chapter  3.  The  model  was  also  used  to  determine  how  long  the  vapor 
delay  at  the  beginning  of  the  postchamber  need  be  to  avoid  transient 
supersaturations.  The  model  also  verifies  that  ice  supersaturation  is 
maintained  throughout  the  postchamber,  so  that  no  crystals  are  lost  by 
subl imation. 

This  model  was  not  run  for  each  experiment  because  it  takes 
approximately  4.5  hours  on  the  Terak.  It  was,  however,  run  for  repre¬ 
sentative  temperatures  and  flow  speeds  used  in  the  studies.  Figure  13 
shows  a  plot  of  the  temperature,  ice  supersaturation,  and  water  super¬ 
saturation  for  the  sample  region  of  a  typical  experiment.  It  is  inter¬ 
esting  to  note  the  behavior  of  the  supersaturation  curves  at  the  begin¬ 
ning  of  the  postchamber.  This  definitely  shows  that  transient  supersa¬ 
turations  are  avoided. 

The  present  model  has  produced  some  results  which  are  somewhat 
different  from  previous  computations.  Substitution  of  the  shear  flow 
profile  for  the  slab  flow  results  in  steady  state  being  attained 
approximately  15X  sooner  but  made  no  change  in  the  steady  state  values. 
This  is  not  surprising  since  as  steady  state  is  achieved  when  and 
go  to  zero,  and  Eqs.  (2)  and  (3)  reduce  to  the  slab  flow  case. 

The  addition  of  total  pressure  and  temperature  dependencies  on  the 
thermal  and  vapor  diffusivities  (Boynton  and  Brattzin,  1929;  List,  1966) 
also  produced  some  changes  in  both  the  transient  time  and  the  steady 
state  values.  The  total  pressure  lowering  results  in  a  significant 


Figure  18.  The  steady  state  temperature  and  supersaturation  with  respect  to  ic 
and  water  for  the  sample  region  through  the  chamber.  The  top  plate 
temperature  is  -10°C  and  the  bottom  plate  -24°C. 
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reduction  in  the  time  required  to  attain  steady  state  conditions  as 
shown  in  Figure  19  but  has  no  significant  effect  on  the  steady  state 
values  themselves.  This  dependency  should  certainly  be  included  for 
laboratories  which  operate  at  atmospheric  pressures  lower  than  the 
standard  1000  millibars  as  well  as  for  any  airborne  measurements. 

The  effects  of  the  temperature  dependency  on  the  coefficients  were 
found  to  be  a  function  of  the  temperature  difference  between  the  plates. 
For  thermal  diffusion  chambers  which  maintain  only  small  temperature 
difference,  as  is  the  case  for  those  used  in  CCN  studies,  the  omission 
of  the  temperature  dependence  is  not  significant.  For  larger  tempera¬ 
ture  differences,  the  effect  becomes  more  apparent  by  making  the  tem¬ 
perature  and  vapor  pressure  profiles  slightly  non-linear.  This  factor 
was  not  considered  by  Katz  and  Mirabel  (1975)  in  their  studies  which 
demonstrated  the  linearity  of  these  profiles.  However,  since  both  pro¬ 
files  are  affected  in  the  same  direction,  there  is  no  significant  dif¬ 
ference  in  the  computed  supersaturation  fields  from  those  using  linear 
temperature  and  vapor  pressure  profiles.  This  effect  in  the  present 
study  lowers  the  temperature  in  the  sample  layer  by  0.2°C  from  the 
linear  temperature  profile  but  produces  no  appreciable  change  in  the 
supersaturation  fields. 

4.3  Experimental  Verification  of  the  Numerical  Models 

There  are  several  experimental  checks  which  can  be  done  to  verify 
the  numerical  models.  First,  since  the  volume  flow  rate  is  measured, 
using  the  chamber  dimensions  and  the  shear  flow  velocity  profile,  the 
maximum  velocity  can  be  computed.  Using  smoke  in  the  sample  layer  and 
positioning  it  midway  between  the  plates,  the  flow  speed  is  measured 
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and  compared  with  that  computed.  A  5  mw  helium-neon  laser  is  used  for 
illumination  of  the  smoke.  Table  2  shows  the  results  of  such  a  compar¬ 
ison  for  the  maximum  velocity  through  the  main  chamber.  The  results 
are  within  limits  of  experimental  errors. 

Since  the  CCN  in  the  atmosphere  are  normally  activated  at  slightly 
above  water  saturation  (Amelin,  1967),  cloud  droplet  should  form  in  the 
chamber  where  the  conditions  are  attained.  Each  of  the  models  can  be 
verified  using  this  method  by  illuminating  the  interior  of  the  chamber 
with  the  laser  beam  and  noting  the  threshold  for  fog  formation.  This 
threshold  has  been  referred  to  as  the  fog  front  (Schaller  and  Fukuta, 
1979).  The  first  model  can  be  verified  by  determining  the  position  of 
the  fog  front  across  the  flow.  The  second  model  can  be  verified  by 
illuminating  the  chamber  along  the  direction  of  the  flow  in  a  region 
where  water  supersaturation  exists  and  determining  the  distance  from 
the  beginning  of  the  chamber  to  serve  within  1.0  cm  of  the  steady  state 
position  indicated  by  the  first  model  and  within  2.0  cm  of  the  tran¬ 
sient  position  indicated  in  the  second  model .  Some  experiments  produced 
"droplets"  below  water  saturation.  These  "droplets"  are  believed  to 
be  haze  particles  instead  of  water  droplets  on  the  nucleated  CCN. 
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TABLE  2 

COMPARISON  OF  THE  COMPUTED  AND  MEASURED  MAXIMUM' 
FLOW  VELOCITY  IN  THE  MAIN  CHAMBER 


Computed  Velocity 

4.2  cm  s'1 

Measured  Velocity 

4.4  cm  s’1 

Relative  Error 

5% 

5.5 

5.6 

2 

6.9 

6.5 

6 

8.6 

8.3 

4 

11.2 

10.3 

9 

13.2 

12.5 

6 

CHAPTER  5 


SAMPLING  PROCEDURE 
5.1  Preparation  of  the  Chamber 

The  preparation  of  the  chamber  for  operation  requires  approxi¬ 
mately  1.5  hours  when  starting  from  room  temperature.  The  procedure 
is  simple  and  does  not  require  the  full  time  attention  of  the  operator. 

Distilled  water  is  used  for  the  ice  surfaces  to  insure  high  accuracy 
values  for  the  boundary  vapor  pressures.  Most  of  these  procedures  were 
discussed  in  the  previous  chapter  but  will  be  briefly  repeated  here  to 
clarify  the  sequence  of  events  in  the  chamber  preparation. 

The  tape  is  applied  to  the  outside  edges  of  the  top  plate  and  the 
predetermined  amount  of  water  is  poured  on  the  bottom  plates  where  ice 
surfaces  are  required.  The  1  cm  wide  copper  plate  is  placed  at  the  be¬ 
ginning  of  the  postchamber.  The  circulating  refrigerated  bath  thermo¬ 
stat  is  set  at  the  desired  top  plate  temperature  and  the  bath  is 
started.  When  the  top  plate  temperature  reaches  -5°C,  the  ice  layer 
is  applied  and  the  tape  removed.  Cooling  of  the  bottom  plate  of  the 
main  chamber  is  begun  with  the  thermoelectric  modules  on  the  cold  side 
drawing  8  amperes  (the  maximum  rating  of  the  modules)  and  those  on  the 
warm  side,  5  amperes  (the  maximum  ouput  of  the  power  supply).  After  the 
water  in  the  main  chamber  freezes,  the  plastic  sheet  to  control  dentrite 
growth  is  inserted.  The  chamber  is  closed,  and  the  top  plate  sealed 
with  the  water  bead,  thermal  insulation  installed  over  the  top  of  the 
chamber.  After  the  bottom  plate  of  the  main  chamber  reaches  the 
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desired  level,  the  thermoelectric  modules  or  the  heating  wire  on  the 
warm  side  of  the  bottom  plate  are  adjusted  to  achieve  the  final  steady 
state  temperature  profile  across  the  plate.  The  thermoelectric  modules 
on  the  cold  side  continue  to  draw  8  amperes  of  electric  current,  i.e. 
continue  at  the  maximum  output.  The  plastic  sheet  is  removed  just 
prior  to  the  beginning  of  sampling. 

5.2  Sampling  Procedures 

The  collection  film  is  placed  in  the  postchamber  and  3  to  6  minutes 
are  allowed  for  the  temperature  to  stabilize.  During  this  time,  the 
water  droplets  condense  on  the  film  with  the  slightly  warmer  top  plate 
being  the  vapor  source.  For  warm  temperatures  (top  plate  less  than 
-10°C),  the  film  is  pretreated  in  the  freezer.  The  air  outlet  device  is 
attached  to  the  postchamber  and  sampling  is  begun  by  adjusting  the  total 
flow  rate  to  the  desired  level  and  then  adjusting  the  sample  flow  rate 
as  required.  The  automatic  temperature  recording  is  started  when  the 
flow  begins.  An  identification  number  anc*  the  date  are  recorded  with 
the  temperature.  The  sampling  continues  until  the  desired  sample  air 
volume  is  processed  in  the  chamber.  The  sample  flow  is  terminated  but 
the  filtered  air  flow  continues  for  15  seconds  to  guarantee  all  the 
sample  air  has  exited  the  chamber.  The  temperature  recording  is  ter¬ 
minated  when  the  flow  is  stopped. 

5.3  Ice  Crystal  Counting 

After  the  sampling  is  completed,  the  cold  cover  with  the  ice  coated 
inner  surface  is  placed  over  the  ice  collection  film  and  both  are  moved 
together  to  the  freezer.  The  film  is  placed  on  a  metal  plate  maintained 
around  -10°C  and  illuminated  for  photographing.  For  crystals  too  small  to 


photograph,  the  wetted  cover  of  1  cm  height  is  placed  over  the  film 
until  the  crystals  grow  to  adequate  size.  Ice  crystals  on  the  film  are 
photographed  with  the  identification  number.  All  the  films  are  photo¬ 
graphed  from  the  same  angle  using  the  same  illumination  source.  After 
photographing,  the  film  is  removed  from  the  freezer  and  dried  at  room 
temperature  for  reuse. 

Three  by  five  inch  prints  are  made  for  each  sampling  case  and  over- 
layed  with  a  grid  showing  the  regions  or  lanes  of  constant  supersatura¬ 
tion  and  temperature.  The  total  number  of  crystals  in  each  lane  is 
counted  and  recorded.  An  average  temperature  is  computed  for  each  of 
the  nine  temperature  locations  for  use  in  the  numerical  models. 

5.4  Sample  Smoke  Generation  and  Collection 

The  artificial  ice  nuclei  used  for  the  experiments  were  generated 
in  the  laboratory.  The  particle  number  conentration  was  determined  by 
the  ultramicroscope  method  (Green  and  Lane,  1964).  It  is  important  to 
determine  the  particle  sizes  and  concentration  (Gerber,  et  al,;  1970). 
The  particle  size  was  estimated  by  measuring  the  particle  fall  velocity 
and  determining  the  particle  size  using  the  Stokes  fall  velocity  equa¬ 
tion  with  the  Cunningham  correction  (Fuchs,  1964).  After  the  number 
concentration  was  determined,  an  appropriate  amount  was  transferred  to 
a  large  box  (approximately  1.5  m  ) to  furnish  a  sample  source  with  a 
reasonable  concentration  and  large  enough  not  to  be  significantly 
affected  by  the  removal  of  sample  for  a  measurement.  Silver  iodide 
smoke  was  generated  by  placing  the  powder  on  a  25  cm  length  of  32  gauge 
platinum  wire  and  heating  the  wire.  The  platinum  wire  requires  42 
watts  for  smoke  generation,  attaining  a  temperature  of  approximately 
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1000°C.  A  photograph  of  the  Agl  smoke  generation  is  presented  in 
Fig.  20(a). 

ON  smoke  was  produced  using  the  device  shown  in  Fig.  20(b).  The 
copper  tube  wound  around  the  soldering  iron  is  heated  and  then  ON 
powder  is  forced  through  the  coil  using  the  syringe.  The  powder  vapor¬ 
izes  as  it  moves  through  the  heated  coil  and  the  smoke  forms  as  the 
vapor  exits  and  mixes  with  outside  air.  The  smoke  is  then  collected 
for  use. 

The  average  diameter  of  the  Agl  particles  was  0-05  pm  3"d  that  of 
DN  0.1  um.  The  size  distribution  are  shown  in  Fig.  21.  The  values  for 
DN  are  from  Vasquez  (1980). 

Natural  nuclei  were  collected  from  air  sampled  from  the  8th  floor 
of  the  Browning  Building,  northeast  corner,  University  of  Utah.  The 
sample  was  brought  into  the  laboratory  through  flexible  steel  plastic 
coated  conduit.  The  kaolinite  sample  was  produced  by  grinding  the 
kaolinite  into  a  fine  powder  and  shaking  it  inside  a  10  liter  plastic 
jug.  Sampling  was  begun  after  a  10  minute  delay  to  allow  large  parti¬ 
cles  to  settle. 


CHAPTER  6 


EXPERIMENTAL  RESULTS 


6.1  Silver  Iodide  (Agl) 

Silver  iodide  was  selected  for  study  since  it  has  been  the  most 
commonly  used  artificial  ice  nuclei  (Flecher,  1962,  Mossop,  1968, 
Cooper,  1970;  Mason,  1971;  Dyer,  1978).  Figure  22  shows  the  results  of 
experiments  under  different  temperatures  and  ranges  of  supersaturation. 
Below  water  saturation,  nucleation  of  water  condensation  does  not  take 
place.  Therefore,  ice  nucleation  that  occurs  under  this  condition  must 
belong  to  the  deposition  nucleation  mechanism.  Since  there  is  little 
change  in  the  activation  mode  of  the  ice  nuclei  across  water  saturation 
up  to  at  least  one  percent  water  supersaturation,  it  appears  that  con¬ 
densation-freezing  is  not  occurring  until  water  supersaturation  condi¬ 
tions  of  at  least  this  value  are  reached. 

Ice  nucleation  is  a  time  dependent  phenomenon  (Warburton  and 
Heffernan,  1964;  Anderson  and  Hallett,  1976)  and  it  is  important  to 
verify  that  the  present  method  of  ice  nucleus  activation  and  counting 
indeed  detects  all  the  active  nuclei  under  the  chosen  condition.  For 
the  condensation-freezing  mechanism  studied  in  these  experiments, 
varying  the  flow  speed  through  the  chamber  can  provide  the  data  needed 
for  this  verification.  If  all  the  time  dependent  nucleus  activations 
are  completed  within  a  period  of  time  shorter  than  available  in  the 
main  chamber,  the  total  number  of  ice  crystals  will  remain  the  same 
regardless  of  the  flow  rate  within  a  reasonable  limit.  Three 
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experiments  were  performed  with  the  same  sample  source  ( Ag I )  under 
identical  chamber  conditions  except  for  varying  flow  speeds.  The  re¬ 
sults  are  presented  in  Fig.  23.  The  average  diameter  of  the  smoke 
particles  is  0.05  um.  For  the  case  where  most  of  the  crystals  were 
collected  on  the  film,  a  near  normal  distribution  was  observed.  The 
other  cases  show  the  same  shape  with  different  half  widths.  If  these 
distributions  are  taken  to  be  normal  and  the  curves  extrapolated  beyond 
the  limits  of  the  collection  film,  then  the  areas  under  the  curves  are 
found  to  agree  within  10%.  This  indicates  that  the  total  number  of 
crystals  nucleated  remain  the  same  for  all  three  cases,  i.e.  there  is 
no  dependence  of  the  total  number  of  ice  nucleations  on  sample  flow 
speed  within  the  limits  of  these  experiments.  It  should  be  observed 
that  the  number  density  of  collected  crystals  falls  to  zero  for  the 
two  slower  flow  speeds  before  the  end  of  the  collection  film.  Thus,  it 
was  confirmed  that  the  necessary  period  of  time  for  all  the  nucleation 
to  take  place  is  shorter  than  that  available  in  the  main  chamber.  For 
deposition  regions  of  these  experiments,  similar  results  were  noted 
although  the  number  of  crystals  was  much  smaller. 

The  experimental  results  shown  in  Fig.  22  clearly  indicate  that 
the  threshold  for  condensation-freezing  is  highly  temperature  dependent 
for  a  given  supersaturation  value.  The  maximum  number  of  nuclei  acti¬ 
vated  is  also  shown  to  depend  on  the  temperature.  For  a  given  tempera¬ 
ture  this  maximum  is  reached  at  some  water  supersaturation  value  but 
becomes  independent  of  it  thereafter.  It  should  be  noted  that  even  for 
the  coldest  temperature,  this  maximum  is  not  reached  until  around  1.5 
percent  water  supersaturation  and  requires  2.5  percent  at  -8.5°C. 

These  high  supersaturation  values  are  rarely  found  in  nature.  Hence, 


of  the  ice  crystal  nunbcr  density  on  the  collection  film  for 
iple  flow  speeds.  The  sample  is  Agl,  the  temperature  is 
10  water  supersaturation  is  2.6%.  The  distance  is  measured 
ling  edge  of  the  collection  film. 


85 


it  is  clear  from  these  data  that  cloud  chambers  which  permit  unrealis¬ 
tically  high  supersaturations  to  occur  do  not  give  a  true  indication  of 
the  activation  characteristics  of  these  nuclei  in  real  clouds  but  in¬ 
stead  are  probably  yielding  the  maximum  values  indicated  in  Figure  22. 

Only  a  very  small  fraction  of  the  particles  are  activated  below 
1.0  percent  water  supersaturation  at  any  of  the  temperatures  studied 
suggesting  that  the  condensation-freezing  mechanism  is  not  yet  present. 
Previous  studies  by  other  investigators  (Isono,  1966;  Schaller  and 
Fukuta,  1979)  have  shown  that  the  condensation-freezing  threshold  is 
lower  than  that  given  here,  the  difference  being  attributed  to  the 
wide  difference  in  instrument  accuracy.  In  addition,  a  very  low  number 
of  nuclei  are  activated  below  water  saturation  in  the  deposition  nucle- 
ation  region.  Measurements  were  also  done  with  the  Agl  smoke  generated 
in  air  filtered  of  CCN  and  IN.  Additionally,  samples  were  used  from 
a  ground  Agl  generator  furnished  by  North  American  Weather  Consultants 
of  Salt  Lake  City.  No  significant  change  from  the  above  results  were 
observed  with  either  sample. 

6.2  1,5  -  Dihydroxynaphthalene  ( DN) 

Experiments  similar  to  those  for  Agl  were  performed  using  1,5  - 
dihydroxynaphthalene  (DN) .  The  results  are  shown  in  Figure  24.  The 
average  size  of  the  smoke  particles  is  0.1  urn.  The  results  are  similar 
to  those  for  Agl  but  differ  in  the  following  respects.  The  water 
supersaturation  threshold  required  for  condensation-freezing  nucleation 
is  lower,  suggesting  that  it  is  more  effective  in  natural  cloud  envi¬ 
ronments.  Again  as  was  the  case  for  Agl,  a  maximum  value  for  the 
active  number  of  nuclei  is  evident  at  high  supersaturations.  However, 
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the  total  number  of  particles  which  are  active  is  greater  by  a  factor 
of  about  three.  For  the  warmer  temperatures  in  Figure  24  the  rate  at 
which  the  number  of  active  DN  particles  changes  with  increased  supersa¬ 
turation  is  slower  than  for  Agl  showing  that  once  the  condensation¬ 
freezing  supersaturation  threshold  is  exceeded,  Agl  reaches  its  maximum 
active  number  more  rapidly  than  does  DN  although  DN  achieves  a  higher 
activation  number. 

6.3  Natural  Nuclei 

Using  this  instrument,  natural  nuclei  have  been  studied  under 
accurately  controlled  temperature  and  supersaturation  without  substrate 
effects.  The  results  are  shown  in  Figure  25.  It  can  be  seen  in  the 
figure  that  the  data  are  considerably  scattered  as  would  be  expected 
from  the  variation  of  natural  aerosol  particle  conditions.  There 
appears  to  be  Tittle  or  no  supersaturation  dependency  for  values  of 
water  supersaturation  between  -1.0  and  1.0  percent  for  temperature 
ranges  from  -17°  to  -22°C.  However,  for  increases  in  water  supersatur¬ 
ation  beyond  1.0  percent,  there  is  indeed  an  overall  increase  in  the 
number  of  active  nuclei.  An  order  of  magnitude  increase  occurs  in  the 
number  as  the  value  of  water  supersaturation  changes  from  1.0  to  2.5 
percent. 

There  were  several  occasions  where  an  interesting  phenomenon  was 
observed.  On  March  26,  1980,  measurements  were  performed  with  outdoor 
air  at  the  northeast  corner  of  the  8th  floor  of  the  Browning  Building, 
University  of  Utah.  The  chamber  temperature  was  -23°C  with  water  super¬ 
saturations  ranging  from  2.9  to  4.0  percent,  the  conditions  where 
active  nuclei  are  normally  observed.  However,  only  three  nuclei  were 
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activated  in  10  liters  of  sample  air.  It  was  later  confirmed  that 
snow  was  falling  at  the  time  of  the  sampling  indicating  that  either 
the  ice  nuclei  in  the  environment  had  previously  been  activated  and 
the  resulting  crystals  excluded  out  of  the  sample  air  or  were  scavenged 
from  the  air  by  the  falling  snowflakes. 

Another  interesting  observation  occurred  on  March  4,  1980.  Out¬ 
door  air  was  being  sampled  at  -20°C  with  water  supersaturation  ranging 
from  1.0  to  3.1  percent.  A  measurement  just  prior  to  a  snowshower 
resulted  in  about  10  active  nuclei  per  liter  in  a  10  liter  air  sample. 
The  next  sampling,  taken  during  the  snowshower,  showed  a  drop  to  less 
than  one  nucleus  per  liter,  i.e.  only  four  nuclei  were  activated  in 
the  10  liter  sample  air.  Since  loading  and  evaporation  of  hydrometeors 
are  the  main  reasons  for  downdraft  development,  it  is  possible  that 
downdraft  air  contains  ice  nuclei  from  evaporated  ice  crystals.  A 
downdraft  prior  to  the  approaching  snowshower  may  have  brought  down 
active  nuclei  which  were  once  involved  in  hydrometeor  formation  process. 
This  sequence  clearly  demonstrates  the  high  variability  of  natural  ice 
nuclei  in  both  time  and  space. 

6.4  Kaolinite 

Since  kaolinite  has  been  studied  by  other  investigators  (Roberts 
and  Hallett,  1967;  Schaller  and  Fukuta,  1979}  and  has  shown  some  abil¬ 
ity  to  act  as  an  ice  nuclei,  several  measurements  were  performed  using 
it.  The  results  are  presented  in  Figure  26.  Although  considerable 
spread  is  seen  in  she  data,  the  same  basic  insensitivity  to  water 
supersaturation  is  evident  as  observed  with  natural  nuclei  until  high 
water  supersaturations  are  reached. 
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Figure  26.  Active  number  of  active  kaolinite  particles  as  a 
function  of  supersaturation  with  respect  to  water 
The  temperature  is  -22°C. 
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CHAPTER  7 

\  CONCLUSIONS 

In  the  attempt  to  develope  a  new  horizontal  gradient,  continuous 
flow,  ice  thermal  diffusion  chamber,  a  number  of  new  schemes  have  been 
devised.  The  subsequent  use  of  the  chamber  for  ice  nucleus  measure¬ 
ments  has  resulted  in  a  number  of  new  findings  both  in  the  technology 
of  instrumentation  for  ice  nucleus  studies  and  in  the  mechanisms  and 
nucleating  behaviors  of  both  natural  and  artificial  ice  nuclei.  The 
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main  contributions  are  listed  below. 

7.1  Instrument  Development 

The  possibility  for  a  high  accuracy  continuous  flow,  horizontal 
gradient,  ice  thermal  diffusion  chamber  was  suggested  by  Fukuta  and 
Saxena  (1979a).  Development  of  the  instrument  presented  in  this  study 
is  a  demonstration  of  the  concept.  Main  features  of  the  developed 
instrument  are  as  follows: 

1.  The  continuous  flow,  horizontal  gradient,  ice  thermal  diffu¬ 
sion  chamber  has  been  constructed  and  successfully  tested.  It 
consists  of  three  sections  which  preprocess  the  sample,  acti¬ 
vate  the  nuclei  under  clearly  defined  conditions  and  collect 
the  nucleated  crystals  for  counting.  In  the  chamber,  the 
sample  air  is  continuously  processed  under  a  range  of  super¬ 
saturations  at  approximately  constant  temperature. 


*1 
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2.  A  scheme  of  sandwiching  the  sample  air  between  layers  of 
filtered  air  and  of  positioning  the  sample  air  in  a  region  of 
clearly  defined  temperature  and  supersaturation  has  been  used. 
By  controlling  the  temperature  and  vapor  pressure  of  the  fil¬ 
tered  air  before  entering  the  chamber,  transie  t  condensation 
in  the  sample  is  avoided  during  preprocessing. 

3.  Heat  pipes  have  been  effectively  utilized  for  temperature 
uniformity  in  the  direction  of  the  flow  in  the  main  chamber. 

4.  For  the  region  where  the  cold  sample  air  enters  the  warmer 
postchamber  after  processing,  a  unique  vapor  diffusion  delay 
scheme  is  applied  to  avoid  transient  supersaturations. 

5.  An  ice  coating  technique  is  developed  to  produce  smooth  ice 
layers  of  uniform  thickness  on  the  inner  surfaces  of  the  top 
and  bottom  plates. 

6.  A  new  ice  crystal  collection  scheme  is  used  which  integrates 
the  ice  crystals  activated  under  the  same  conditions.  Uncer¬ 
tainties  associated  with  the  optical  methods  used  in  other 
diffusion  chambers  are  eliminated. 

7.  Use  of  filtered  air  to  separate  the  sample  from  the  walls 
eliminates  any  wall  effect  on  the  sample.  Moreover,  it  keeps 
the  sample  air  away  from  the  convectively  unstable  region 
adjacent  to  the  walls. 

8.  The  flow  control  scheme  employed  permits  accurate  monitoring 
and  control  while  introducing  the  sample  without  any  contamin¬ 
ation  from  pumps  or  flowneters. 


I 


9.  The  condensation-freezing  and  deposition  mechanisms  of  ice 
nucleation  have  been  effectively  isolated  while  the  contact 
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freezing  mechanism  has  been  excluded.  This  feature  permits 
the  identification  of  the  mechanism  involved  in  ice  nuclei 
activation. 

10.  A  scheme  for  high  accuracy  computerized  real-time  temperature 
measurement  has  been  developed.  This  permits  precise  deter¬ 
mination  of  the  interior  temperature  and  supersaturation 
fields. 

11.  Complete  numerical  models  have  been  constructed  which  furnish 
detailed  temperature  and  vapor  pressure  fields.  The  models 
not  only  compute  steady  state  values  but  verify  that  transient 
supersaturations  and  wall  effects  are  avoided. 

7.2  New  Findings  with  Ice  Nuclei  Measurements 

The  accurate  temperature  and  supersaturation  control  available 
with  the  present  instrument  has  enabled  the  findings  of  new  as  well  as 
more  detailed  behaviors  of  both  artificial  and  natural  ice  nuclei.  The 
main  findings  are  listed  below. 

1.  The  rate  of  deposition  nucleation  is  much  lower  than  that  of 
condensation-freezing  for  all  nuclei  in  all  temperature  ranges 
studied. 

2.  For  Agl  of  average  particle  size  0.05  pm,  only  a  very  small 
fraction  (less  than  0.001%)  are  found  to  be  active  below  1.0 
percent  water  supersaturation  and  for  DN  of  average  particle 
size  0.1  urn,  a  higher  number  (as  high  as  0.02%)  are  active. 

3.  For  temperatures  warmer  than  -10°C,  the  condensation-freezing 
mechanism  is  not  effective  for  Agl.  This  finding  can  have  a 
significant  implication  on  the  use  of  Agl  in  weather  modifica¬ 


tion. 
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4.  The  total  number  of  active  ice  nuclei  at  high  supersaturations 
is  about  three  times  larger  for  DN  than  that  for  Agl. 

5.  Once  the  condensation-freezing  thresholds  are  reached  for  Agl 
and  DN,  further  increases  in  supersaturation  makes  Agl  to 
reach  its  maximum  active  number  more  quickly  than  DN. 

6.  Natural  and  kaolinite  nuclei  show  little  supersaturation  de¬ 
pendency  below  1.0  percent  water  supersaturation. 

7.  Natural  and  kaolinite  nuclei  show  an  increase  in  activity  by 
an  order  of  magnitude  as  the  water  supersaturation  increases 
from  1.0  to  2.5  percent.  Instruments  which  allow  such  high 
supersaturation  values  to  occur  will  therefore  yield  high 
active  nuclei  counts  and  the  results  should  be  considered 
unrel iable. 

8.  For  natural  nuclei,  large  variances  were  observed  with  high 
counts  in  cloud  downdrafts  and  low  counts  during  precipitation. 

7.3  Possible  Future  Applications  of  the  New  Chamber 

The  new  features  of  the  present  chamber,  i.e.  the  accurate  and 
simultaneous  display  of  a  range  of  supersaturation  with  continuous 
sampling,  suggest  application  of  the  chamber  in  a  number  of  areas. 

1.  The  present  chamber  may  be  used  as  a  standard  for  detecting 
condensation-freezing  and  deposition  nucleations  and  other 
chambers  and  methods  of  ice  nucleus  detection  may  be  calibrated 
against  it. 

2.  Further  detailed  studies  of  condensation-freezing  and  deposi¬ 
tion  nucleations  may  be  carried  out  for  different  chemicals, 
particles  of  different  sizes,  effects  of  additives  on  ice 


nucleus  smoke  and  nucleation  time  lag. 

The  present  chamber  may  be  used  to  obtain  real  time  data  of 
natural  nuclei. 

Contact  freezing  mechanism  may  be  isolated  out  of  mixed 
nucleation  processes  using  the  present  chamber. 


APPENDIX  A 


NUMERICAL  MOOEL  TO  DETERMINE  WALL  EFFECTS  AND 
THE  STEADY  STATE  SUPERSATURATION  FIELDS 

This  numerical  model  computes  the  steady  state  temperature  and 
vapor  pressure  values  inside  a  rectangular  thermal  diffusion  chamber 
using  the  same  approach  presented  by  Tomlinson  and  Fukuta  (1979).  The 
model  provides  for  temperature  and  vapor  pressure  gradients  and  there¬ 
fore  the  supersaturation  profile  inside  the  chamber  utilizing  thermal 
gradients  on  both  the  top  and  bottom  plates  although  in  this  study, 
only  the  bottom  plate  maintains  a  gradient.  This  generality  allows  the 
model  to  be  used  for  other  thermal  diffusion  chambers.  Also,  the  model 
is  not  limited  to  subfreezing  temperatures  and  can  be  used  for  any  CCN 
spectrometer  of  rectangular  geometry  such  as  that  of  Fukuta  and  Saxena 
(1979a  and  b) . 

The  model  assumes  linear  temperature  gradients  in  the  vertical 
direction  on  the  side  walls  and  throughout  the  chamber  interior  as 
justified  by  Tomlinson  and  Fukuta  (1979).  On  the  other  hand,  any  tem¬ 
perature  profile  can  be  imposed  on  the  top  and  bottom  plates.  In  this 
study,  the  top  plate  is  kept  isothermal.  The  temperature  is  recorded 
during  the  experiments  at  three  locations  across  the  bottom  plate  of 
the  main  chamber.  The  temperature  profile  for  the  bottom  plate  is 
taken  to  be  a  quadratic  function  of  the  distance  horizontally  across 
the  flow  direction  which  satisfies  these  three  temperature  points. 

Vapor  pressure  is  assumed  to  be  ice  saturated  on  all  the  boundaries 
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(water  saturation  for  temperatures  above  0°C)  and  is  computed  using  a 
sixth  degree  polynomial  (Rasmussen,  1978). 

The  equations  used  are  taken  from  Carslaw  and  Jaeger  (1957).  In 
order  to  use  these  equations,  the  boundaries  have  to  be  expressed  by  a 
sine  series.  A  sufficient  number  of  terms  are  used  to  express  exact 
values  at  the  boundary  grid  points.  The  following  series  were  used: 


n  /  \  r  ■  nrrx 

Vx)  '  L  «„  s'"  — 

n=u 


p,(x)  =  I  b  sin  T  ,  and 
L  n=0  n  B 


n-rrz 


PRU)  ■  l  c„  sin  -j 
n=0 


where  an>  b  ,  and  cn  are  the  amplitudes  of  wave  number  n  in  each  series 
(Spiegel,  1968).  The  lower  left  corner  of  the  chamber  is  the  origin  for 
the  coordinate  system  as  shown  in  Figure  Al.  The  x  direction  is  hori¬ 
zontal  and  has  a  maximum  at  z  =  A  on  the  right  boundary.  The  z  direc¬ 
tion  is  vertical  and  has  a  maximum  at  z  =  B  on  the  top  boundary. 

As  discussed  in  Chapter  4,  the  vapor  pressure  can  be  expressed  by 
the  following  sum  as  shown  in  Figure  A2: 

P(x,z)  *  P]  +  P2(x,z)  +  P3(x,z)  +  P4(x,z)  +  Pg(x,z) 

The  equations  used  to  compute  each  pressure  component  are  as  follows: 


ii  * — i  8  H  8  =1  t> 


The  coordinate  system  used  in  the  model  and  the  boundary 
vapor  pressures. 
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The  values  can  be  computed  on  any  grid  desired.  The  grid  used  in 
this  study  was  16  X  11  which  gives  a  1  cm  resolution  horizontally  and 
0.06  cm  vertically  for  the  chamber  dimensions.  After  the  temperature 
and  vapor  pressure  fields  are  computed,  supersaturation  values  with 
respect  to  both  water  and  ice  are  determined  (water  only  for  tempera¬ 
tures  above  0°C).  These  values  are  compared  with  the  supersaturation 
values  for  the  same  boundary  conditions  in  the  infinitely  wide  parallel 
plate  case  and  the  error  fields  are  estimated  accordingly.  It  is 
these  error  fields  which  are  used  to  determine  the  extent  of  the  wall 
effects.  The  supersaturation  values  for  the  interior  of  the  chamber 
are  used  in  the  studies. 

The  error  associated  with  the  wall  effect  is  presented  in  Figure 
A2.  Several  cases  were  computed  for  various  chamber  heights.  It 
should  be  noted  that  the  wall  effect  increases  significantly  as  the 
chamber  height  increases. 
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Figure  A2.  Error  In  the  water  supersaturation  field  due  to  the 
wall  effect  for  a  rectangular  ice  thermal  diffusion 
chamber.  Various  heights  in  centimeters  are  plotted. 


APPENDIX  B 


SCHEMES  USED  IN  THE  NUMERICAL  INTEGRATION 


As  presented  in  Chapter  4,  the  equations  describing  the  tempera 
ture  and  vapor  pressure  profiles  between  the  plates  of  the  thermal 
diffusion  chamber  with  flow  are 


(Bl) 

(B2) 


The  velocity  profile  for  shear  flow  between  parallel  plates  is 

v<2>  *  -  2?  57  2  <H-2> 

3  P 

Here  V(z)  is  a  known  function  of  z  for  any  given  pressure  r— .  Both  < 

°y 

and  D  are  functions  of  total  pressure  and  temperature  (List,  1966). 
For  a  given  altitude,  i.e.  atmospheric  pressure,  <  and  D  become  func¬ 
tions  of  temperature  only.  Hence,  for  every  time  step  <  and  D  are 
determined  at  each  interior  grid  point  after  the  temperature  at  that 
grid  point  is  computed. 

3T  3  P 

The  time  derivatives  ^  and  are  approximated  with  the  second 
order  centered  approximations  as  follows: 
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where  n  is  the  index  for  time,  j  the  index  for  the  y  direction,  k  the 
index  for  the  z  direction,  and  At  the  time  increment.  The  second  de¬ 


rivatives  are  approximated  by  the  second  order  scheme  as  follows: 
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where  Az  is  the  vertical  grid  spacing.  The  first  derivatives  in  space 
are  approximated  using  second  order  centered  differencing  except  for 
the  right  (downstream)  boundary.  The  scheme  is  as  follow: 
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where  Ay  is  the  horizontal  grid  spacing  along  the  direction  of  the  flow. 
On  the  right  boundary,  the  following  upstream  scheme  is  used: 
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P.  -  P 

V(z)  |y  =  V ( 2 )  (B12) 

For  the  beginning  of  the  postchamber  where  the  vapor  pressure  is 
delayed,  the  new  lower  boundary  condition  produced  some  oscillations 
in  the  solution  field.  The  use  of  the  upstream  differencing  scheme  in 
the  advection  term  resulted  in  a  stable  solution  field.  Since  this  is 
only  a  first  order  scheme,  the  horizontal  grid  spacing  is  reduced  by  a 
factor  of  ten  in  order  to  maintain  the  overall  accuracy. 

The  vertical  grid  spacing  is  z  =  0.06  cm  and  the  horizontal  grid 
spacing  is  y  *  0.5  cm  except  as  mentioned  above  at  the  beginning  of  the 
postchamber  where  it  is  decreased  to  0.05  cm.  These  grid  spacings 
proved  adequate  to  resolve  the  conditions  inside  the  chamber.  This 
was  verified  by  decreasing  the  grid  spacing  by  a  factor  of  two  and 
noting  no  change  in  the  solution  field.  The  time  increment  was  selected 
to  satisfy  the  following  stability  requirements: 

maximum  ]V[At  <Q^  (B13) 
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Az 
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Although  each  of  these  critria  is  sufficient  only  in  particular  more 
idealized  equations,  satisfaction  of  the  most  stringent  of  (B13,  B14, 
B15)  produced  stable  solutions. 
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ABSTRACT 


It  Is  well  recognized  that  conditions  In  the  environment  surround¬ 
ing  Ice  nucleu  (IN)  particles  must  be  accurately  controlled  In  order  to 
understand  their  nucleatlon  behavior.  Concerning  the  condensation¬ 
freezing  and  deposition  mechanisms  of  Ice  nucleatlon,  the  main  factors 
representing  the  environment  are  supersaturation  and  supercooling 
(temperature).  Starting  front  the  concept  of  the  wedge-shaped  Ice  ther¬ 
mal  diffusion  chamber  of  Schaller  and  Fukuta,  and  that  of  the  Fukuta- 
Saxena  cloud  condensation  nucleus  (CCN)  spectrometer ,  and  by  addition¬ 
ally  considering  operational  efficiency  and  accuracy,  a  continuous-flow, 
horizontal  gradient,  ice  thermal  diffusion  chamber  has  been  developed. 
The  chamber  consists  of  three  parts,  l.e.,  preprocessing,  main  activa¬ 
tion  and  Ice  crystal  settling  sections.  A  common  top  plate  Is  used  for 
all  three  sections  and  is  held  Isothermal  using  a  circulating  bath  while 
a  temperature  gradient  Is  maintained  across  the  bottom  plate  of  the  main 
section  by  applying  thermoelectric  modules.  In  the  preprocessing  and 
main  sections,  both  plate's  are  covered  with  ice.  A  newly  developed 
method  Is  applied  to  smoothly  coat  Ice  on  the  plates.  This  configura¬ 
tion  of  the  main  section  results  in  a  nearly  constant  temperature  horl- 
zontally  in  the  direction  across  the  sample  flow,  and  produces  a  range 
of  supersaturations.  Heat  pipes  are  utilized  on  the  sides  of  the 
bottom  plate  to  insure  temperature  uniformity  along  the  direction  of 
the  flow.  The  sample  volume  is  sandwiched  In  the  region  of  maximum 
jupersaturation  between  layers  of  filtered  and  predried  air  such  that, 
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as  the  sample  enters  the  main  section,  a  nearly  constant  supersaturation 
Is  achieved  vertically  through  the  sample  without  transient  supersatura¬ 
tions.  The  preprocessing  section  Is  held  Isothermal  at  the  top  plate 
temperature.  The  design  and  flow  of  the  main  section  permit  the  nucle¬ 
ated  Ice  crystals  to  be  carried  Into  the  Ice  crystal  settling  section 
without  loss.  The  bottom  plate  In  the  settling  section  Is  maintained  at 
a  temperature  slightly  lower  than  that  of  the  top  plate.  Formation  of 
transient  supersaturations  at  the  entrance  of  the  settling  section  has 
been  avoided  by  delaying  vapor  diffusion  while  allowing  thermal  diffu¬ 
sion  to  proceed.  The  problem  of  transient  supersaturation  development 
has  been  examined  for  shear  (Polseullle)  flow  cases  using  a  numerical 
method  Instead  of  the  commonly  used  slab  flow  assumption.  Wall  effects 
have  also  been  computed  for  the  supersaturation  distribution.  The  sta¬ 
bility  of  the  sample  flow  throughout  the  entire  chamber  has  been  con- 
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firmed  with  smoke  tests'.  A  new  method  for  Ice  crystal  detection  has 
been  developed  and  is  employed  in  the  Ice  crystal  settling  section. 

Mylar  copy  film  (carbon  paper)  holding  condensed  water  droplets  Is 
placed  in  the  settling  section.  Ice  crystals  nucleated  In  the  main 
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section  fall  on  the  film  and  grow  to  visible  sizes  in  the  presence  of 
the  droplets.  The  positions  ofthe  Ice  crystals  across  the  flow  direc¬ 
tion  give  the  supersaturations  and  temperatures  at  which  they  nucleated. 
The  optimum  design  of  the  chamber  In  terms  of  the  height,  width,  length 
and  flow  speed  has  been  determined. 

Two  numerical  models  have  been  developed  to  compute  the  temperature 
and  supersaturation  fields  within  the  chamber.  The  first  computes  the 
steady  state  conditions  In  a  vertical  plane  perpendicular  to  the  flow 
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and  Includes  the  wall  effects.  The  second  computes  the  steady  state 

conditions  in  a  vertical  plane  parallel  to  the  flow  using  the  shear  flow 
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velocity  profile  and  the  total  pressure-  and  temperature-dependent  co¬ 
efficients  for  vapor  and  thermal  diffusion. 

The  Increased  accuracy  of  the  developed  chamber  has  resulted  In 
several  significant  new  findings.  For  all  the  temperature  and  supersa¬ 
turation  ranges  studied,  the  rate  of  deposition  nucleatlon  1$  much 
lower  than  that  for  condensatlon-freeting.  For  silver  Iodide  (Agl) 
particles  of  average  diameter  0.05  ym,  only  a  small  fraction  are  active 
below  1%  water  supersaturation,  while  for  1,5  -  dlhydroxynaphthalene 
(ON)  of  average  diameter  0.1  ym,  a  much  larger  fraction  are  active. 

The  total  number  of  active  nuclei  for  these  samples  was  approximately 
three  times  greater  for  DN  than  for  Agl  at  high  supersaturations.  For 
temperatures  warmer  than  -10°C,  the  condensation-freezing  mechanism  Is 
not  effective  for  Agl.  Natural  nuclei  and  kaolinlte  show  little  super- 
saturation  dependence  below  IX  water  supersaturation  and  as  much  as  an 
order  of  magnitude  increase  as  water  supersaturatjon  Increases  from  1% 
to  2.5%.  Additionally  for  natural  nuclei,  large  variances  wer'e  observed' 
with  high  counts  in  cloud  downdrafts  and  low  counts  during  precipitation. 
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